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Abstract 



The thesis is considering aspects of SU{2) Yang-Mills thermodynamics in its deconfining high- 
temperature phase. 

We calculate the two-point correlation function of the energy density of the photon in a ther- 
malized gas, at first in the conventional U{1) gauge theory, followed by a calculation, where the 
photon is identified with the massless gauge mode in deconfining SU{2) Yang-Mills thermody- 
namics. Apart from the fact, that this calculation is interesting from a technical point of view, 
we can consider several aspects of phenomenological relevance. Since we interpret the two-point 
correlator of energy density as a measure for the energy transfer, and thus for the electromag- 
netic interaction of microscopic objects, such as atoms immersed into a photon gas, we are able to 
give an explanation for the unexpected stability of cold, innergalactic clouds consisting of atomic 
hydrogen. 

Subsequently, we evaluate the spatial string tension in deconfining SU{2) Yang-Mills thermo- 
dynamics, which can be regarded as measure for the magnetic flux through the area enclosed by 
the associated Wilson loop. On the level of on-shell polarization eff'ects for the massless mode we 
observe a perimeter-law, and wc speculate that the lattice-obtained area-law is induced by off-shell 
contributions to the polarization tensor. Moreover, we discuss an interesting two-loop result for 
the pressure which seems to be associated with the presence of screened magnetic monopoles being 
responsible for an area-law. 



Zusammenfassung 

In dieser Diplomarbeit werden Aspekte der SU{2) Yang-Mills Thermodynamik in ihrer dekonfi- 
nierten Hochtemperatur-Phase betrachtet. 

Wir berechnen die Zwei-Punkt Korrelationsfunktion der Energiedichte des Photons in einem 
thermalisiertem Gas, zunachst in der konventionellen [/(l)-Eichtheorie. Danach betrachten wir 
den Fall, dass das Photon mit dem masselosen Eichfeld in einer dekonfinierten SU{2) Yang-Mills 
Theorie identifiziert wird. Neben der Tatsache, dass diese Rechnung vom technischen Standpunkt 
aus von Interesse ist, konnen wir mehrere Aspekte phanomenologischer Natur betrachen. Sofern 
wir den Zwei-Punkt Korrelator der Energiedichte als ein Mafi fiir den Energietransfer, und damit 
auch fiir die elektromagnetische Wechselwirkung mikroskopischer Objekte, wie in ein Photongas 
eingebettete Atome, ansehen, sind wir in der Lage, eine Erklarung fiir die unerwartete Stabilitat 
kalter, innergalaktischer Wolken, bestehend aus atomarem Wasserstoff, zu geben. 

Anschliefiend berechnen wir die raumliche Stringspannung in dekonfinierter SU{2) Yang-Mills 
Thermodynamik. Diese Grosse kann als Mai3 fiir den magnetischen FluB durch die vom ent- 
sprechenden Wegner- Wilson Loop umschlossene Flache interpretiert werden. Auf dem Niveau von 
on-shell Polarisationseffekten fiir die masselose Mode beobachten wir ein Perimeter-Gesetz, und 
wir vermuten, dass das aus Gittersimulationen bestatigte Flachen-Gesctz auf off-shcU-Bcitrage des 
Polarisationstensors zuriickzufiihren ist. Weiterhin diskutieren wir ein interessantes Zwei-Schleifen 
Resultat fiir den Druck, das mit der Prasenz abgeschirmter magnetischer Monopole in Zusammen- 
hang zu stehen scheint, welche wiederum verantwortlich fiir ein Flachen-Gesetz waren. 
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Chapter 1 

Introduction 



The importance of gauge theories in theoretical physics is evident, since they provide a deep, 
mathematical framework for the formulation of quantum field theories. Fundamental interactions 
are described by gauge theories. Interactions, as they occur in the standard model of particle 
physics (SM), between fermionic fields are mediated via the exchange of bosonic gauge fields. The 
interaction is completely determined through gauge symmetry. In particular, non-abelian gauge 
groups proved to be an important step in describing Nature on the level of elementary particles. 

In the SM gauge theories are considered in the framework of perturbation theory. Many results 
derived on this basis are verified up to enormous precision. For examples we refer to the anomalous 
magnetic moment of the electron in Quantum Electrodynamics (QED) or the asymptotic freedom 
in Quantum Chromodynamics (QCD). Nevertheless, it is not the case that perturbative methods 
leave no open questions. At first we have to mention that not even the most elementary processes 
of QED are mathematically described in a sound way. Perturbative solutions are asymptotic series, 
the philosophy is an expansion in the small gauge coupling. Thus the perturbative concepts are 
only viable for small couplings. In case of a strongly coupled system, the philosophy of perturba- 
tion theory collapses. For instance, the dependence of the gauge coupling on the energy scale is 
characterized by the /3- function. In QCD, for increasing energy, the coupling decreases logarithmi- 
cally and thus the calculation of the /3-function in terms of a perturbation series seems to be well 
justified for large momentum transfer, see also [T]. This phenomenon is called asymptotic freedom 
and was discovered in 1973 by Gross, Wilczek and Politzer 0. As for the thermalized situation 
perturbation theory leads to infrared instabilities, resulting from weakly screened soft magnetic 
modes (5- This manifests itself in a gluon plasma at high temperature, and is in contradiction to 
naive expectation, since one would expect perturbative methods to yield reasonable results, be- 
cause the gauge coupling approaches zero for temperature T exceeding the Yang-Mills scale many 
times over. However, a perturbative expansion of the thermodynamical pressure in a Yang-Mills 
gas merely is applicable up to order g^hig. 

Beyond this order we have to regard the fact, that in pure non-abelian gauge theories exist 
infrared stabilizing, topologically nontrivial field configurations in the euclidean formulation, while 
in perturbation theory the expansion is about the trivial vacuum. 

In addition, there are a number of experimental and observational results in elementary particle 
physics and cosmology, where the SM either yields no explanation or even leads to contradiction 
with the experiments and observation. For instance, there is the fact that the SM does not provide 
for any theory about the phenomenon of dark matter and dark energy, which corresponds to about 
95% of the total energy density content of the Universe. Furthermore the direct experimental 
detection of the Higgs particle is still open. This leaves room for speculations concerning the 
electroweak symmetry breaking of the SM. Another unexplained observation is the missing power 
in the CMB temperature-temperature correlation at large angular scales [U [6] • Finally we refer 
to the existence of a large, cold, old, dilute, innergalactic cloud (GSII130-03-69) in-between spiral 
arms of the Milky Way [71 [S]. The astrophysical origin of this and similar structures is mysterious, 
its high content of atomic hydrogen in view of its unexpectedly large age clashes with the SM, see 

^We recall that already in 1969 I. Khriplovich 2 and in 1972 G. 't Hooft ^ had performed according calculations. 



1 



SecEH 

In order to overcome these problems an investigation of (non-abelian) gauge theories involving 
nonperturbative methods seems to be essential. In [S] an analytic and nonperturbative approach 
to SU (2) and SU (3) Yang-Mills thermodynamics was developed. The construction of this effective 
theory is close in spirit to the Ginzburg-Landau-Abrikosov theory, a macroscopic, phenomenono- 
logical description of superconductivity in metals, where the microscopic processes leading to a 
nontrivial ground state are summarized by scalar fields and pure-gauge configurations. In the 
deconfining phase of SU{2)/SU{3) Yang-Mills thermodynamics [l^ the thermal ground state is 
determined by an adjoint scalar field and a pure-gauge field configuration a'^f , solving the Yang- 
Mills equation in the presence of the background </). This macroscopic field originates from spatial 
coarse-graining over microscopic processes. The field 4> is compose d of noninteracting calorons (and 
anticalorons) of topological charge one and trivial holonomy. This composite adjoint scalar field 
governs the dynamics of the spatially averaged, BPS saturated, topological nontrivial and noninter- 
acting part of the thermal ground state, it is quantum mechanically and statistically inert, and thus 
it indeed acts as a background for the topological trivial sector of the theory, which is responsible 
for caloron interactions. The sector of trivial topology is minimally coupled to the coarse-grained 
caloron configurations and manifested as the already mentioned pure-gauge configuration a^J^. The 
macroscopic adjoint scalar acts as a Higgs field, and it breaks the SU{2) gauge symmetry down 
to C/(l). Thereby, on tree-level, two out of three gauge bosons acquire a temperature-dependent 
mass, while the third gauge mode remains massless. Another remarkable property of the effective 
theory is its finiteness in the ultraviolet and infrared. An infrared cutoff results from caloron in- 
duced gauge boson masses, an ultraviolet cutoff is provided by the compositeness scale |0| which 
constrains momenta in loop calculations in the only physical gauge: unitary-Coulomb. Thus no 
renormalization procedure to absorb divergencies is necessary. The evolution of the effective gauge 
coupling e is derived from the requirement of thermodynamical self-consistency. That is, Legendre 
transformations following from the fundamental partition function, need to be kept intact after 
spatial coarse-graining. In the deconfining phase of SU{2) Yang-Mills thermodynamics the effective 
coupling has a plateau value e{T) — ^/Stt, and diverges logarithmically at the critical temperature. 
This is the phase boundary to the preconfining (magnetic) phase. In this thesis we consider the 
high-temperature deconfining (electric) phase of SU{2) Yang-Mills thermodynamics only. 

The thesis is organized as follows: Chapter two gives an introduction into the concepts of 
Yang-Mills theories. Nonperturbative, topological field configurations in euclidean spacetime are 
presented, especially instanton, or at finite temperature, caloron configurations. The next section 
contains a review of the developments in nonperturbative, analytic SU (2) Yang-Mills thermody- 
namics. The generation of the nonperturbative ground state and the construction of the effective 
theory is presented, followed by a brief overview of the other two phases of the theory. Chapter 
three contains the calculation of the two-point correlation function of energy density, at first in 
the conventional U{1) electrodynamics, followed by a calculation in deconfining SU{2) Yang-Mills 
thermodynamics. We compare our results for both cases and observe a sizeable suppression of 
the SU (2) correlator at low temperatures and large distances. Since we interpret the two-point 
correlation of energy d ensity as a measure for the energy transfer and thus for the interaction 
between single atoms, immersed in a photon gas, we give a possible explanation concerning the 
unexpected stability of cold clouds consisting of atomic hydrogen located in-between spiral arms 
of our galaxy. Chapter four is dedicated to a quantity called spatial string tension. This object 
is defined as the logarithm of the spatial Wilson loop and is a measure for the dual magnetic flux 
through the area enclosed by the loop. In lattice gauge theory the high-temperature properties 
of the string tension were calculated. On the level of on-shell 'photon' polarization effects the 
generally accepted lattice results are not reproduced, but we are able to offer an explanation why 
our approach is insufRcient at this level. Finally, in Chapter five we give a summary of the thesis 
and an outlook on further research. In the appendix, we recall some general concepts of quantum 
field theory at finite temperature and state the Feynman rules in the real-time formulation. In 
addition, we give some technical details concerning various calculations in this work. 
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Chapter 2 

Yang-Mills Thermodynamics 



2.1 Lie groups and Lie algebras 

The following section is dedicated to the mathematical requirements for the formulation formulation 
of non-abelian gauge theories. We give a short introduction into the theory of Lie groups and Lie 
algebras. More material concerning the theory of Lie groups can be found in [TT | fT ^ IT5]. 

Lie group: A Lie group is a differentiable manifold G endowed with an associative, differentiable, 
invertible map 

GxG^G (51,52)^51-52 V5i,52 e G, (2.1) 

such that G obeys the group properties. The existence of a neutral element e and an inverse 
element 5^^ is provided. 

As examples for Lie groups we regard the group of regular invertible linear transformations of 
a vector space V , Gl{V). Of particular interest in theoretical physics are the Lie groups Gl{N,l$.) 
and Gl{N, C) and some of their subgroups. It is a crucial fact that every subgroup of a Lie group is 
a Lie group itself. If V is given as the A^-dimensional euclidcan space R^, the general linear group 
Gl{N,M.) consists of all real N x A^-matrices with dctg 7^ 0. Gl{N,M.) is a manifold of dimension 

If V is given as the complex vector space , the Lie group Gl {N, C) is the set of all nonsingular 
A^ X A^-matrices with complex entries. As subgroup of physical relevance we should mention the 
A^^ — 1-dimensional special unitary group SU{N), consisting of all unitary matrices M G U{N) 
with unit determinant. 

The center of a group G is the subgroup containing all those elements x in G, that commute 
with all other elements, that is xg = gx for all g in G. 

Lie algebra: A Lie algebra is a vector space together with a bilinear operator [•, •] 

VxV^V {x,y)^[x,y] yx,yeV, (2.2) 

called the Lie bracket, which is anticommutative, 

[x,y]=-[y,x] Vx,yeV, (2.3) 
and satisfies the Jacobi identity 

[x,[y,z]] + [z,[x,y]] + [y,[z,x]]^0 yx,y,z£V. (2.4) 

If we constrain ourselves to Matrix groups, the Lie bracket is identical with the commutator of two 
matrices. For a Lie group Gl{V) the associated Lie algebra is denoted as 2l{V). 

As an example for a Lie algebra we point out that the tangent space at an arbitrary element 
g of a Lie group G always possesses the structure of a Lie algebra. We assume that the group 
elements 5 G G are parametrized by a set of continuous parameters a" (a = 1, 2, ... , A") in a way 
that a = corresponds to the identity, i.e. 5(a° = 0) = e. The set of parameters represent a 
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differentiable manifold. Since the group elements depend analytically on the parameters a", any 
element g £ G (at least in the vicinity of the identity element) can be Taylor expanded about the 
unit element e of the group. It holds 



where t"^ 



dg{a) 
da"- 



(2.5) 



Q = 



The objects are referred to as infinitesimal generators of the Lie group. The generators of a Lie 
group always establish a Lie algebra, where the dimension of the vector space V (and with this in 
mind the dimension of the Lie algebra) corresponds to the number of generators K . In order to 
illustrate this fact we consider the concrete case of a group manifold, parametrized by coordinates 
x^, x^, . . . , x'^ . A basis for the tangent space is spanned by the directional derivatives di = A 
tangent vector v^di acts on a differentiable scalar function / defined on G as 



v(/) 



(2.6) 



Of fundamental importance is the tangent space T^^G at the unit element of G and the Lie bracket 
of two tangent vector fields x and y, acting on a differentiable scalar function / on G, defined as 



[x,y\ = x{y{f)) - v{x{J)) , 



(2.7) 



building up a Lie algebra. 

Since the generators constitute a basis, the Lie bracket of two generators is proportional to 
another basis vector 



(2.8) 



where the real and antisymmetric factors f"''"^ are called structure constants of the algebra. For 
the Lie algebra su(2) a possible basis is given by the Pauli matrices cr^''^-^ multiplied with the 
imaginary unit «, these three skew hermitian matrices are generators of 5u(2). The commutator 
relation (|2.8p is in this particular case 



T' y 



ahc 



(2.9) 



where e"^"^ denotes the Levi-Civita tensor. Analogously su(3) is spanned by the set of Gell-Mann 
matrices. 

A Lie algebra is connected to its associated Lie group via the exponential map, defined as 
the power series exp(^) = Yl^=o matrix A is an element of the algebra, then exp(^) 

is an element of the group. Thus two conditions on the matrices are implied: A ought to be 
anti- hermitian {A^ = —A) and traceless {trA — 0). 

The rank of a Lie algebra is then defined as the dimension of the Cartan subalgebra, which is 
given as the maximal algebra of commuting generators. 

Linear representation: A linear representation of a Lie group G on a vector space V is a ho- 
momorphism which connects every element x G G with a linear transformation R(x) acting on V . 
The dimension of the representation R{x) is equal to the dimension of the vector space V. If a 
basis of V is assigned, R{x) can be expressed as a matrix. 

A representation is called (ir)reducible if the representation can(not) be transformed into the 
blockdiagonal form 



R{x) 



R\x) A{x) 

i?2(x) 



(2.10) 



If A(x) vanishes for any x E G, R{x) is called fully reducible. 

A linear representation of the Lie algebra is formed by matrices which satisfy the commutation 
relation (|2.8p of the Lie algebra. Since the elements of a matrix group are linear transformations, 
the group elements themselves generate the fundamental representation of the group. 

The r generators of a Lie group provide a representation of dimension r. A basis of this adjoint 
representation is then provided through the structure constants, Rabtc — ifabc- The structure 
constants obey the Jacobi identity and the commutation relation as required. 
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2.2 SU(N) Yang-Mills theory: Foundations and topological 
aspects 



In this section we consider the formulation of non-abeHan gauge theories: The Lagrangian density 
and the equation of motion for the gauge fields is derived. We give a review of solitonic field 
configurations solving the Yang-Mills equation. As a first example for an object of topological 
origin we motivate the emergence of magnetic monopoles in thermalized, pure Yang-Mills theory 
according to the construction by 't Hooft and Polyakov. Subsequently, the main focus lies on a 
solution to the equation of motion in euclidean field theory, which is stabilized by topology, called 
instanton. This is followed by a presentation of its finite temperature generalization, the caloron. 
An introduction into foundations of non-abelian gauge theories and topological aspects of quantum 
field theory can for instance be found in [HI [HI [HI [T7] . The topic of instanton physics is covered 
in [m [19] on a fundamental level. The generalization of instantons to finite temperature field 
theory is reviewed in [20l [21] . In this thesis we use the convention that Greek indices correspond 
to four-dimensional Lorentz indices, while the Latin indices (from the middle of the alphabet) 
denote spatial indices. Lie algebra indices are denoted with Latin indices from the beginning of the 
alphabet. The Einstein summation convention is implied if not stated otherwise. In Minkowskian 
spacetime we have to distinguish between upper and lower indices. Only contractions between 
covariant and contravariant indices are possible. In Euclidean spacetime, a distinction between 
upper and lower indices is not necessary. 



2.2.1 The Yang-Mills Equation 

In the following section we contemplate a four-dimcnsional Minkowskian spacetime. In 1954 Yang 
and Mills |22j constructed a Lagrangian density, which is invariant under non-abelian gauge trans- 
formations. This Lagrangian describes a nontrivial interacting field theory. While in the original 
work the gauge symmetry was taken to be SU (2) of isospin rotations, soon a generalization to arbi- 
trary non-abelian gauge groups took place. The perturbative analysis of these theories was a merit 
of Feynman, Faddeev, Popov, DeWitt, and many other authors in the 60's and 70's, when Feyn- 
man rules for their perturbative quantization were developed. The tremendous physical relevance 
of non-abelian gauge theories as a mathematical framework for the description of fundamental 
interactions is generally acknowledged. 

We designate an SU{N) Yang-Mills theory as being pure, if only gauge fields and no matter 
fields are involved. We demand from the Lagrangian of an SU (N) Yang-Mills theory gauge invari- 
ance under any local SU{N) transformation U{x). In order to construct gauge-invariant kinetic 
quantities from the gauge fields in the Lagrangian, all partial derivatives, acting on non-gauge 
fields in a non-pure Yang-Mills theory, are replaced by covariant derivatives, defined as 

D^ = d,,-igA,,, (2.11) 

where is a gauge field and g denotes the gauge coupling of the fundamental theory. This gauge 
field is an element of the Lie algebra, and therefore an expansion exists in terms of the generators 
A° of the algebra 

A,^Al^. (2.12) 

The generators are normalized as 

trA^A** = (2.13) 

For determining the effect of a covariant derivative on a non-gauge field 0, we have to consider the 
representation of the field. If 4> transforms under the fundamental representation we obtain 

D^c^^d^cj^-igAl^c^, (2.14) 
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which is different from the case when (f> is in adjoint representation. Then we have 

iD^^)^d^(t>~ig[A^,^]. (2.f5) 

Since we require that the covariant derivative is constructed in a manner, so that the derivative of 
the field obeys the same transformation law as the field itself, the gauge field transforms under a 
gauge transformation U {x) £ SU {N) according to 

A^{x) -> U{x)A^{x)U\x) + -U{x)d^{x)U\x) . (2.16) 

In the language of differential geometry the gauge field can be regarded as a connection on an 
analytic manifold, the gauge group. The field strength tensor is associated with the curvature on 
the manifold and is defined through the Lie bracket of covariant derivatives: 

F^, = ^[D^,D,]^d^A,-d,A^~ig[Af,,A,], (2.17) 

and in component notation 

= d^K - + ar'^AlAl . (2.18) 

In non-abelian gauge theories the field strength tensor is not gauge-invariant, but transforms 
homogeneously 

F^,,{x) ~* U{x)F^,{x)U\x) . (2.19) 

Therefore, the quantity ivF^^F^^'^ is a gauge invariant Lorentz scalar, and the Lagrangian density 
is given as 

£--itrF^,F^'' = ~ii^;,F'^''^ (2.20) 

where the normalization of the generators (I2.13P was used. In the case of Abelian gauge groups the 
commutator vanishes and therefore no self-interaction of the gauge bosons occurs. In contrast, for 
non-abelian gauge groups self- interactions governed by three- and four-boson- vertices are present, 
due to the fact that the Lagrangian contains additional cubic and quartic terms in the gauge 
fields. Although a mass term for gauge fields is for reasons of gauge invariance not admissible, the 
gauge bosons can acquire a mass by dynamical symmetry breaking via the Higgs mechanism, if 
the accordingly charged scalar sector is added to the Lagrangian in Eq. (|2.20p . 

The Euler-Lagrange equations are derived from minimizing the action associated with the 
Yang-Mills Lagrangian ((T^ . We get 

D^.F'^" = , (2.24) 

or in component notation 

a^F^"" + gf'^bc^i.bptj.vc ^ {^D^pt^^Y ^ . (2.25) 

According to Eq. (|2.25p , the fact that a pure SU (N) Yang-Mills theory is interacting is reflected 
in the nonlinearity of the field strength tensor and the term proportional to g. 



^To construct a theory of gauge fields interacting with fermions we have to add the Dirac Lagrangian to the 
gauge-field Lagrangian H2.20| l: 

C = -jFI'^^F>""^ +iiPlpil>-miiip. (2.21) 

The Euler-Lagrange equations are obtained by varying the resulting action and we obtain 

(j0-m)i/) = O, (2.22) 

as the Dirac equation for the fermion field, and 

{D>'F^,r = -9ji, (2.23) 

as equation of motion for the gauge field, where = 'tp^ji^t'^ iIj is the current density of the fermion field. 
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2.2.2 The 't Hooft-Polyakov monopole 

An abelian J7(l)-gauge theory of electromagnetism containing magnetic monopoles in addition to 
electric charges is more symmetric than the standard Maxweh theory. Despite that, there is no ne- 
cessity for the existence of magnetic monopoles. However, if the gauge theory of electromagnetism 
would be described by an underlying non-abelian gauge group (like SU(2)) subject to dynamical 
symmetry breaking, the field equations would yield solutions with magnetic monopole constituents. 
The existence of magnetic monopoles at finite temperature is enabled by the nontrivial homotopy 
group Tr2{SU{2)/U{l)) = Z [531 [511 [5S]. The theoretical possibility of magnetic monopoles in 
gauge theories was elucidated the first time when 't Hooft [26] and Polyakov [27] considered an 
S'C/(2)-Yang-Mills theory with an isovector adjoint Higgs field (/>"'. The Lagrangian density is then 
given as 

1 1 77?^ 

^ = -iK'^F""' + ^{D^r){D^r) - —rr - x{rrf ■ (2.26) 

If the parameter m? is chosen negative, the vacuum expectation value of the Higgs field minimizes 
the action at = =; p. 

In 1974 't Hooft derived non-singular solutions with finite energy to the field equations resulting 
from the Lagrangian (|2.26p by considering static configurations. Setting Aq = 0, the SU{2) gauge 
field and the Higgs field has the following asymptotic behavior (|x| — > oo): 

At = -ef^ (2.27) 

r = F^^. (2.28) 

It is worth mentioning that in these configurations space and isospace indices are mixed. For a 
field configuration of the form (j2.27p the name 'hedgehog' was introduced by Polyakov. 

The SU (2) gauge symmetry is broken due to the nonvanishing vacuum expectation value of 
the Higgs field. The field strength is determined by the S'[/(2)-gauge invariant 't Hooft-tensor 
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•^M- = - zrn^e'^'^r{D^<f>'')iD,^n , (2.29) 



where |0| = ^J4>°-(j)°-. Assuming that the gauge field and the scalar field point into the three direction 
of the Lie algebra, we insert 

4=^^^, ^^ = F^0. (2.30) 

Fixing the scalar field (/) in isospace in such a way, the 't Hooft tensor reduces to the usual abelian 

U{\) electromagnetic field tensor for the gauge field A^^. Inserting the hedgehog- Ansatz (|2.27|) into 
Eq. (|2.29p leads to an asymptotic static field strength 

■^y = "^j^^??'^" ■ (2-31) 
This corresponds to a radial magnetic field 

Ba = ^. (2.32) 

The magnetic flux follows from 
f Air 

^ = i daaBa = — , (2.33) 

and the magnetic charge as 

= ^=e-^ (2.34) 
47r 
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In the BPS-limit (to ^ 0, A ^ 0) the monopole mass is evaluated as Mm = ^Mw, where Mw 
denotes the vector boson mass, given as Mw = e\4>\- We define the current 

= = lef^-^^d^F^x = ~^e''''^^e'''''d,rdj''dx$'' , (2.35) 

2 2e 

where = 

The magnetic charge can then be written as the integral over a sphere at infinity: 

— / (fxK° = - 

Att / Sire 



q^ = — I d^xK" = -— (p da' £'^'=e"''=<^"5,/(9j/ . (2.36) 



52 

This equation leads to two important statements. The magnetic charge depends only on the 
asymptotic behavior of the fields and the magnetic current merely has contributions from the Higgs 
field. The current is identically conserved, its divergence vanishes. Notice that this current 
conservation is a consequence of topology, and does not follow from a continuous symmetry of the 
Lagrangian. 

2.2.3 Instanton solution 

In order to construct a four-dimensional well localized solution to the gauge field equation in the 
case of a pure SU{2) Yang-Mills theory, we have to consider this theory in euclidean spacetime. 
Thus we perform a Wick-rotation in the zero-coordinate, t — ir, where r € M. The Yang-Mills 
action is then given as H 

Jd^x tiF^.F^, - ^ Jd^x tr(E2 + B^) , (2.37) 

with the color electric and magnetic fields E'* and B'*. A solution to the equation of motion is 
determined by detecting the minima of the action. The existence of the configuration is guaranteed 
by 7r3(5'C/(2)) = Z. The SU{N) Yang-Mills vacuum is degenerate, that is there exist apart from 
the vacuum = infinitely many configurations of the kind 

A^^iUix)d^uHx), (2.38) 

which differ from the trivial vacuum only by an SU{N) gauge transformation U{x), and thereby 
have zero field strength and energy momentum tensor. These so-called pure gauges fall into distinct 
topological classes, which are characterized by an integer winding number 

nw = J d^xe''''tr{U^d^U){U^djU){U^dkU) . (2.39) 

There exists no continuous transformation between two pure gauges with different winding number, 
so a transition can only be realized via a tunneling process with probability amplitude e~'^. Such 
a tunneling solution is called instanton. 

We define the dual field strength tensor as 

F^:y — ~£fii^KxFfj,i/ ■ (2.40) 

In the dual field strength the electric fields are replaced by the magnetic fields, and vice versa. 
Now we are able to write the Yang-Mills equation in the form of a Bogomolnyi decomposition 

S = J d'^xtrF^.F^, 



^ J d^xtv \±F^,F^, + -{F^, T F^,){F^, T F^,) ] . (2.41) 



^In the remainder of this section we absorb the gauge coupling into the gauge field. This is a common convention 
for considering nonperturbative aspects. 
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The first summand in tlie Yang-Mills action (I2.4ip is the topological charge also denoted as Pon- 
tryagin index 

Q=T^ fd^^i^ F^^uF^u I d^x tr . (2.42) 



16^2 J ^A'-'^- 

The integrand is called Pontryagin density and is equal to the total divergence of the Chern-Sinions 
current Kn. It holds thalU 



d^K^ - ^F;,F^,, (2.48) 



1 

32^2-''''"'"" 

with 



(2.49) 



The winding number defined in ()2.39p is connected to the charge of the Chern-Simons current by 
virtue of 



nw = j d'xK° = j^J d'xe^.k [A-^d^Al + i/'^^M^^A^ ) . (2.50) 



Notice that the Chern-Simons current is not a gauge invariant quantity, in contrast to the Pon- 
tryagin density. Provided that the integrand is nonsingular and rapidly decreasing in the limit of 
spatial infinity, Gauss' theorem yields again the topological charge 

dr-^ f d^xKn+ I dr I da^ Ki 



dr . 

d^x Ko+ / d^x Kq 



= nw{T = oo) — nw{T = — oo) . (2-51) 

We recognize that a tunneling field configuration from one vacuum to another topologically dif- 
ferent is associated with a Pontryagin index Q 7^ 0. Thus the second term in the Bogomolnyi 



''Consider for instance the case of SU{2) as gauge group |16l 117) . The Chern-Simons current can be written 



i^M = T^^'"'«^ [Kd.Al + -e^^-A^A^^A^j , (2.43) 

where e"*"^ denotes the Levi-Civita tensor. We obtain tliis form of by employing the commutator relation 
^■^j -^j = ie"^'^^ and the normalization condition of the generators of the Lie-algebra tr -^-^ = ^'^"''i where 
^a,b,c denote the Pauli matrices. Thus can be written as trace 



= ^^M-^Atr i^-A^d^Ax - -A^A^Axj . (2.44) 
With respect to the cyclicity of the trace the divergence of this quantity is then given as 

d^K^ = K^ = -^£„,«A tr (l-{dtiA^){d^Ax) - iid^A,,)A^Ax) . (2.45) 



47r2 

Now we calculate tv F^yF^i, = \e^vK.\Pti.vF^x i with _Fp^ = flpA^ — d^A^ — i[A^,Ay] . We obtain 

tr F^„F^„ = 2e^^^x tr {d/iAv){d^Ax)-2ie^^KX tr {d^Ai,)A^Ax-2ie^^KX tr A^A^(aKAA)-2e^^KA tr A^Ai,AkAx ■ 

(2.46) 

Due to the cyclicity property of the trace the second and the third terms are equal and the last term vanishes for 
reasons of symmetry. Thus in the case of SU (2) 

^ K.K. (2.47) 



" 327r' 
is proven. 
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decomposition (|2.4ip is a square, and the euchdean action is bounded from below (Bogomolnyi 
bound). The minimum in a given topological sector is reached for an selfdual or a anti-selfdual 
field configuration 

= ±P;. ■ (2.52) 

A selfdual field corresponds to an instanton, an anti-selfdual field to an antiinstanton. The mini- 
mum of the action is then given as 5 = We mention that the nonperturbative character of 
instantons is reflected in the essential zero at g = in the transition amplitude, thus pcrturbative 
expansions completely disregard instanton solutions. 

The selfduality-equation ()2.52p constitutes a differential equation of first order, in contrast to 
the second order like Euler-Lagrange equation (|2.24|) . The Bianchi identity, 

D^F^x + DxF^^ + D^F^x = 0, (2.53) 

provides that (anti-)selfdual field configurations satisfy the equations of motion automaticalljQ. 
(Anti-)sclfdual configuration have a vanishing energy-momentum tensor. This is also called BPS- 
saturation[29j. 

In order to construct an explicit solution to ()2.52|) we have to demand that the field strength 
vanishes at infinity. This condition is realized, if the gauge field approaches a pure gauge solution 
on the boundary of euclidean spacetime. 

A gauge field which satisfies the mentioned requirements, should show an asymptotic behavior 

like 

(N-«^), (2.54) 
where \x\ denotes ^x^x^j, and the definition of the 't Hooft symbols is given as 

'nl.u = ^atiu + 5a^i.6vA " 5avS^,A , (2.55) 
V% = ^a^iu - 5a^5vi + 5av5^A ■ (2.56) 



Eq. (|2.54p is generalized to the case of finite |x| as 

Al{x)^2r^a^..xJ-^, (2.57) 
where / fulfills the boundary conditions 

lim /(x^) -> 1 lim /(a;^) ^ const ■ x^ . (2.58) 

a:^ — *oo — >0 

The request for selfduality determines a differential equation for /, which subject to Eq. (|2.58p is 
solved by 

/=^^. (2.59) 

The constant of integration p is called the instanton radius. The translational invariance ensures 
that the center of the instanton can be shifted to an arbitrary point . 



1 
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Summarizing we obtain the Belavin-Polyakov-Schwarz-Tyupkin instanton [30j 



{x - z)y 



with the corresponding field strength 



[{X - zf + p2] 



(2.60) 
(2.61) 

^a^^v ■ This field configuration is well localized 



The antiinstanton is obtained by substituting T^a^;^ , ,iaf^ 
in spacetime, since the field strength decreases like x~'^. 

The topological charge of the instanton is derived via evaluation of Eq. (|2.5ip with the asymp- 
totic Chern-Simons current 



hm A„, = — . 



\x\- 



27r2 x^ 



(2.62) 



The BPST-instanton carries the topological charge Q = 1 (respectively Q = — 1 for the BPST- 
antiinstanton). 

The BPST-instanton (|2.60p has eight free parameters (also denoted as collective coordinates): 
four due to the translation of the instanton center, the instanton radius, associated with dilata- 
tions and three coordinates in color space. The action of the BPST-instanton is independent of 
these parameters, so each of them gives rise to a zero mode, i.e. they are moduli of the instanton 
solution. 

The BPST-instanton so far was considered in regular gauge. This implies that the configuration 
is well behaved except at infinity where a singularity emerges. However, this singularity can be 
shifted from infinity to the center of the instanton by application of a suitable gauge transformation. 
In this singular gauge the gauge field and the field strength transform to 



2?7a^,.(a; - xq)^ 



= -Vaf^i^diy In 



P 



{x — Xo)^ [{x — z)2 + p2J 
„2 



{x - xoY 
jx - xo)^(x - Xo)p ^ 
{x - Xo)"^ 



Xo) 



(2.63) 
(2.64) 



In the gauge potential there now occurs a singularity at the instanton center, but nevertheless 
the field strength remains finite. The fact that in singular gauge the action density as well as the 
topological charge Q is localized at the instanton center enables us to construct field configurations 
with topological charge higher than unity, so-called multi-instanton solutions. This construction 
can be realized by generalization of the singular gauge BPST-instanton to the 't Hooft-ansatz for 
the gauge field 

A^ix) ^ ~f]a^,dAryTiix) , (2.65) 

where 11(2:) denotes a pre-potential. The required selfduality of the field strength yields a 
Laplace-like equation for the pre-potential 



d^d^Il{x) 



= 0. (2.66) 



n{x) 

A general solution with topological charge Q = K evaluates to 

K 
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i.e. it describes instantons with centers at Zi and radii pi. (The muhi-instanton solution with 
Q — K is constructed by substituting fja^iv through r]a^v in the 't Hooft-ansatz (|2.65p . A muhi- 
instanton with Pontryagin index Q = K can be considered as a superposition of K single instantons, 
thus parameters generate the moduli space. But we should mention that the pre-potential 
p.67p does not represent the most general solution, since all K instantons have the same orientation 
in the color space. The most general construction of a multi-instanton was obtained in 1978 by 
Atiyah, Drinfeld, Hitchin and Manin [5T] . 



2.2.4 Caloron solution 

Quantum field theory at finite temperature ([X]) can be formulated in euclidcan signature. Thereby 
temperature T is identified with the euclidean time r compactified on a circle with perimeter 
/3 = . The fields have to show /3-periodicity, at least up to a gauge transformation 

A^(/3, x) - (7(0, x)A^(0, x)C/t (0, x) + Jc/(0, x)5,,(0, x);7^(0, x) . (2.68) 

The intention is then to construct an instanton solution which obeys this condition. The /9-periodic 
field configuration is called caloron, the /3-periodic anti-selfdual solution anticaloron [32] . 

The initial point concerning this construction is the 't Hooft-ansatz (|2.65p with a pre-potential 
as defined in (|2.66p . We obtain the caloron centered at (0, 0) associated pre-potential by evaluating 
the infinite sum over instanton pre-potentials centered at (n/3, 0),n G Z, and equal radii pi. In 
1978 Harrington and Shepherd calculated this quantity 

^ r? -Kf? sinh ( 2f 1 

-oo^^~^^^'^> cosh (^r) ~ COS (^1 

where r — |x|. Integrating the associated Chern-Simons current over a three-sphere (with center at 
the singularity), we obtain the topological charge of the caloron (anticaloron) as plus one (minus 
one). The action of the Harrington- Shepherd caloron calculates as equal to the instanton action 

o 2 

An on the classical level temperature- independent value. 

In analogy to the construction of multi-instanton configurations in singular gauge a caloron of 
topological charge higher then unity is built as a superposition of the pre-potentials of the invoked 
single calorons. Such a multi-caloron solution possesses among the radius more dimensionful 
moduli. However, they do not yield any contribution to the thermodynamics, as we explain below. 

Calorons are categorized in field configurations with trivial and nontrivial holonomy. The 
essential quantity is the value of the Polyakov loop at spatial infinity, which is a topological 
invariant. The Polyakov loop is defined as a time- like Wegner- Wilson loop at finite temperature 

ig f drj44{r.x) 

P{x)^Ve , (2.70) 

where V is the path-ordering operator. Calorons with a Polyakov loop lying in the center of the Lie 
group are said to have trivial holonomy, otherwise they have nontrivial holonomy. In particular, for 
the interesting case of SU{2) calorons with P(|x| oo) = ±1 corresponds to trivial holonomy. For 
the Harrington-Shepherd caloron we obtain for the Polyakov loop P — I (in the limit \x — > oo|). 

In the early 80's Nahm [23^ pointed out the existence of caloron solutions with nontrivial 
holonomy, in 1998 an explicit construction was accomplished indepedently by Lee and Lu ^24, and 
Kraan and van Baal |25_ . Lee and Lu parametrized 

A^(t,|x|-^(X3) = -*|A3, (2.71) 

where A3 is the skew-Hermitian generator and u denotes the holonomy, which is constrained to 
< It < The solution of the ADHMN-equations [23l [3T] yield a selfdual field configuration 
consisting of BPS-magnetic (anti-)monopole constituents. This is the nontrivial holonomy caloron. 
The masses of the constituents are given as 

TOl = inu , 7712 = Ivr - 7?j . (2.72) 
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The one-loop effective action S^a (or the one-loop quantum weight) of a trivial-holonomy caloron 
was computed by Gross, Pisarski and Yaffe [2U], while the calculation in the nontrivial-holonomy 
case was performed by Diakonov, Gromov and Slizovskiy in 2004 [33j . Summarizing those results, 
we can conclude that the kind of interaction between BPS monopoles and antimonopoles in the 
caloron is depending on its holonomy. Small holonomy leads to an attractive potential, large 
holonomy to a repulsive potential. Calorons with nontrivial holonomy are no stable objects, since 
they either decay into a into a monopole-antimonopole pair or they annihilate each other and, 
then subsequently the caloron collapses back onto trivial holonomy. In contrast to the emergence 
of small holonomy calorons, the genesis of a large holonomy caloron (^u ~ is strongly Boltzmann- 

suppressed subject to e-'^(™i+™^) ~ e-^'''. As a consequence the ground-state physics is essentially 
governed by small holonomy calorons. 



2.3 The deconfining phase of SU(2) Yang-Mills thermody- 
namics 

In this section we give an outline of [5] where fundamental concepts of SU(2)-Yang-Mills thermo- 
dynamics are developed. An investigation of the microscopic dynamics of individual calorons is a 
sophisticated task, and for the derivation of macroscopic quantities as pressure and energy density, 
not convenient if not impossible. But as long as we are focussing on the thermodynamics a spatial 
average via coarse-graining over the physics resulting from the topological nontrivial sector is a 
feasible method, which yields exact results. 



2.3.1 Construction of the composite adjoint Higgs field 

The spatial average over calorons leads us to the concept of a macroscopic thermal ground state, 
which is characterized by a macroscopic field (f> and a pure-gauge configuration a^*. The adjoint 
scalar field (j) is generated by calorons with topological charge one and trivial holonomy. We require 
for the following properties: 



1. Due to spatial isotropy must be a Lorentz scalar. 



2. In a pure Yang-Mills theory all local fields (as elements of the Lie algebra) transform under 
the adjoint representation of the gauge group. So as composite of this fields transform 
under the adjoint representation as well. 

3. Since describes a homogeneous ground state of a thermal system its modulus is independent 
of space and time. This modulus is governed by a dynamically generated Yang-Mills scale A 
and temperature. 

4. The color orientation of (p, also referred to (/)'s phase, is r-dependent. Since (j) is composed of 
calorons the phase of (p is periodic in euclidean time. The classical caloron action S = ^p- is 
independent on temperature, so (/)'s phase is not explicitly time independent. 

Thus the field can be written as a product of modulus and phase 



In the following we sketch the evaluation of 0's modulus and phase and present the results. 
(f> 's phase: In [3 [31] it is shown that the phase of the field can be evaluated as follows: 



tr / d^x / dpA'^F^,(r,0){(T,0),(r,x)}^^,,,(r,x){(r,x),(T,0)}, (2.74) 
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Here the field strength belongs to a Harrington-Shepherd caloron with caloron scale parameter p, 
and 



{(T,0),(r,|x|)} 
{(T,|x|),(r,0)} 



-tr<^2 ^ 



= V cxp 



V exp 



i J dZf^Af_,{z) 

(r,x) 



(2.75) 
(2.76) 

(2.77) 



The Wilson lines are calculated along straight lines. The ~-sign declares that ^ as well as the 
trace in equation ()2.74p solve the same homogeneous evolution equation in r 



V 



101 



= 0, 



(2.78) 



where V represents a linear second-order differential equation and turns out to be uniquely deter- 
mined as 



(2.79) 



Eq. (|2.74p is to evaluate on both caloron and anticaloron field configuration and the results have 
to be added up. Indeed, it transforms like an adjoint scalar. We underline that the construction of 
(^'s phase is well justified and reasonable, since every generalization leads to explicit temperature 
dependencies or additional scales. 

1. Due to the selfduality of calorons and anticalorons every local definition of 0's phase yields 
zero. 

2. Inclusion of higher n-point functions (n > 3) would force us to introduce additional factors 

order to get a dimcnsionless contribution. However this would lead to an explicit 
/3-dependence of the action on the classical level. 

3. Inclusion of nontrivial-holonomy calorons is not admissible, because of the divergence of the 
one-loop effective action in the thermodynamical limit 20]. 

4. A shift of the initial point of the Wilson line from (r, 0) to an arbitrary point (r, z) would 
introduce a finite mass scale |z|^^ which is not existent on the classical level. 

5. The calculation of the Wilson lines needs to be performed along straight lines, because spatial 
curvature would correspond to an additional mass-scale, which would have to be compensated 
by a factor of /3. A non-flat measure of p has the same effect. 

Another condition to the macroscopic field is a consequence of BPS-saturation. Since (f) is con- 
structed from non-interacting BPS-saturated calorons, the energy-momentum tensor vanishes iden- 
tically for the composite object. Thus (p solves a first-order differential equation, the EPS-equation 



(2.80) 



where y(i/2) a 'square-root' of a the potential V^cj)) = tr{V^^/'^^)^V'^'^/^K This BPS-equation is 
determined from the evolution equation for the phase up to global gauge rotations and a choice of 
winding sense. In order to evaluate the phase we have to find a traceless, hermitian configuration 
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which solves the second-order equation of motion (|2.78p for the phase as well as the first-order 
BPS-equation (|2.80p . In [HlIMj a solution in the (1, 2)-plane of the Lie algebra is derived as 

^ = CiAiCOs(^^(T-ro)) +C2A2sin(^^(r-ro)) (2.81) 
/ 27ri 

= CAiCxp ( =F — A3(t-to) 

with undetermined integration constants. Since the phase is r-dependent, we refer to this gauge as 
winding gauge. This solution corresponds to circular polarisation in the X1-X2 plane of color space. 
(j) winds along an on the group manifold of SU{2). Notice that a global gauge ambiguity 
associated with the choice of the polarization plane of this oscillation remains. 

(f) 's modulus: Initially, we assume the existence of an externally given, constant Yang-Mills scale 
A. A is a renormalization group invariant. This is only realized by an effect resulting from trivial- 
topology fluctuations called dimensional transmutation. If we require that the phase Eq. (|2.8ip is 
reproduced, ]/(^/^) ought to be a linear function of (f>. The temperature-dependence of (p cannot 
be explicit, since the average over the caloron-anticaloron system is /3- independent. But V((j)) may 
depend on temperature through the periodicity of 4>. Furthermore, the potential should be an 
analytic function of (p. The conditions lead us to a BPS-equation for (j) 

d,cl> = ±iA^X^^. (2.82) 
Inserting the expression for the phase, we obtain for ^'s modulus O [31] 



H(/3,A) = V^ = V^ = A, (2.83) 

where A denotes the dimensionless temperature A = 27rTA~"'^. Thus the field (j> vanishes propor- 
tional to T~2 with increasing temperature. An important fact is that the modulus corresponds 
to the maximal resolving power allowed in the effective theory. 

Since 0's phase and modulus are determined we are able to calculate the potential V{(j)) from 
the BPS-equation Eq. (jTSUl) as 

F((/)) = A^tr0~2 = AttTA^ . (2.84) 

The euclidean Lagrangian governing the dynamics of the Higgs field is then given as 

£0 = tr ((3,0)2 + AV"') • (2.85) 

We consider now the statistical and quantum mechanical inertness of (/> [F. A field is called inert, 
if the ratio of the squared masses of fluctuations Sip, defined as ^5— an the squared compositeness 
scale 101, or respectively the temperature, are both larger than unity. We obtain 

^ = 12A3 and = 487r^ (2.86) 

In the deconfining phase of SU{2) Yang-Mills thermodynamics it holds A > Ac = 13.87, and the 
criterions for inertness are satisfied. We conclude that quantum fluctuations of the field do not 
exist, since their masses would exceed the compositeness scale many times over. 



2.3.2 The effective theory 

At first we recall the validity of a unique decomposition for every SU{2) gauge field of the fun- 
damental theory. Each gauge field Af^ decays into a topologically nontrivial, BPS-saturated part 
A^°P represented by calorons and a topologically trivial remainder a^. Until now, we have confined 
ourselves to the topologically nontrivial sector, described by the adjoint scalar field p after spatial 
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average. In order to include caloron interactions we have to consider the topological trivial (i.e. 
Q — 0) sector. The topologically trivial fluctuations are responsible for a change of short-lived 
holonomy of calorons and for interactions between magnetic monopole constituents of nontrivial 
holonomy. In this situation the inert scalar field 4> acts as a background after coarse-graining. 
The influence of the configurations a^u is taken into account via minimal coupling 

d^,(j)-^ D^,(f)^d^,(l)-ie[af,,(p], (2.87) 

where e denotes the gauge coupling in the effective theory (which is not equal to the coupling g of 
the fundamental theory). 

We are interested in an effective theory, after the application of spatial coarse-graining. Integra- 
tion of topological defects into an inert field (j) and the renormalizability of Yang-Mills theory 
assures form-invariance of the effective Lagrangian including the topologically trivial sector under 
the applied spatial coarse-graining. The resulting gauge-invariant, effective Yang-Mills action de- 
termining the dynamics of Q = 0-fluctuations and the background field cj) in the deconfining 
phase is then given as 

5[a^] ^ JdrJ d?xtr Qg^.G^. + [D^cjyf + , (2.88) 







where Gf^i, — G°^— = d^au — d^a^ — ie[a^^ a^]. The Euler-Lagrange equation concerning the field 



flp follows as 

D^G^,^ie[4>,D,<t)\. (2.89) 

For reasons of rotational invariance this equation of motion requires a pure-gauge solution a^^ , 
implying Di^cj) = [1]. It holds that 

flS; = l(a^J7)17t ^ with n = e±*lf . (2.90) 

At this point, the ground-state is completely determined by a^'* and cj). The energy-momentum 
tensor associated with this ground-state is then 

= 47rA3r,5^, , (2.91) 

implying a ground-state equation connecting pressure and energy density: 

pa^ ^ -At:K^T = -/^ . (2.92) 

Since the ground state contains not only contributions of the BPS-saturated fields, but as well from 
topological configurations, P^'' and p^^ attain finite values after spatial-coarse graining. Thereby 
the Yang-Mills scale becomes gravitationally measurable. 

In the following we discuss propagating fluctuations 5a^ in the effective theory. It is convenient 
to transform from winding gauge to unitary gauge via the transformation 

[/^g-ifA2^_ (2.93) 

In unitary gauge the pure-gauge configuration vanishes, while the adjoint scalar field is rotated 
into the three-direction of the algebra: 



0, 0=|0|A3. (2.94) 



This gauge transformation switches the value of the Polyakov loop from Vla^f — (5^4^ A3] = — 1 
to Vla^f = 0] = 1. Therefore we can conclude that the theory possesses two equivalent ground 
states. The fact that the ground state is degenerate with respect to a global electric Z2{SU2) center 
symmetry designates the electric phase to be deconfining. (j) can be beheld as an adjoint Higgs field, 
it breaks the SU{2) gauge symmetry dynamically down to the Abelian subgroup U (1). Due to the 
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(1 2) 

adjoint Higgs mechanism two out of the three gauge modes Saji ' acquire a temperature-dependent 



mass 



A3 A 

m = 2e{TM = 2e{T)\l — = 2e(r)-^ , (2.95) 

while the third gauge mode remains massless. Here e(T) is again the temperature-dependent 
effective gauge couphng, which wiU be investigated in detail below. The massive bosons are referred 
to tree-level heavy (TLH) and are denoted as V'^ and V~ , the massless boson as tree-level massless 
(TLM) and denoted as 7 . 

Finally, we review the emergence of compositeness constraints on the maximal off-shellness of 
quantum fluctuations. The existence of a maximal resolution |(/)| resulting from spatial coarse- 
graining introduces limits for the propagating momenta. Momenta exceeding the compositeness 
scale 101 would lead to a higher resolution and are thus forbidden. We state the constraints: 

1. The maximal off-shellness of a gauge mode is restricted to 

|p2 _ < |0|2 ^ (2.96) 
where p and m denote four-momentum and mass of the quantum fluctuation. 

2. The four- momentum flux in a four- vertex is constraint subject to 

\{p + kf\<\cf>\\ (2.97) 

In this case two of the four legs have to form a loop, otherwise a differentiation between s-, 
t- and u-channel scattering is mandatory. 

Since the phase space for fluctuations is tremendously limited, the massive excitations are 
very weakly interacting thermal quasi-particles. That implies the negligibility of higher-loop con- 
tributions, thus (at large temperatures) the order of magnitude from loop-order to loop-order is 
about 10"'^ — lO""*. It is worth mentioning that the generated masses yield infrared cutoffs, and 
the arising of infrared divergencies in loop expansions is avoided automatically. Furthermore, as 
a side-effect the compositeness constraints avoid the emergence of ultraviolet divergencies, so no 
standard renormalization procedure is necessary in the effective theory. 

The evolution equation for the effective gauge coupling e(T) is derived from the demand for 
thermodynamical self-consistency. Thermodynamical quantities are connected via Legendre trans- 
formations in the fundamental theory. Self-consistency is provided, if a reformulation of the theory 
into an effective Lagrangian does not alter the Legendre transformation interrelating the thermo- 
dynamical quantities. Consider the Legendre transformation which maps the total pressure onto 
the total energy density: 

dP 

P = T—-P. (2.98) 

The requirement of self-consistency implies that the pressure shows only explicit temperature 
dependence. Since the massive gauge modes obtain a temperature dependent mass, the condition 
can be formulated in terms of a dimensionless mass a = (3m = /3e|0| = 27reA~'^/^ as 

daP = Q. (2.99) 

In [5] the one-loop expression for the pressure was calculated as 

= -A^ [2P(0) + 6P(2a)] + 2a) , (2.100) 

where the function P{a) is defined as 



P(a)= I dxx'ln 1 _ Q-Vx^+a- _ (2.101) 
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Employing the self-consistency condition Eg . (|2.99p we get the evolution equation for A (a): 

2A\^a D{2a) 
WTT2prb(2a) 

where the function D{a) is defined as 



Dig) ^ I dx , — — . (2.103) 



Inversion of the differential equation yields an evolution equation for the effective gauge coupling 
e(A). 

2.3.3 Phase transitions: The preconfining and the confining phase 

In this section we summarize the essential features of the other phases in the theory of SU{2) 
Yang-Mills thermodynamics. A comprehensive representation of the physics in these phases of the 
theory can be found in ,9, ,3^ [371 [3E| ■ 

In the deconfining phase of SU{2) Yang-Mills thermodynamics e(A) is mostly constant at a 
plateau value at e = -v/Stt, except in the vicinity of the critical temperature at Ac = 13.87, where 
the gauge coupling diverges logarithmically 

e(A) - -log(A- Ac). (2.104) 

At this point a phase transition to the preconfining (magnetic) phase occurs. 

The preconfining phase: The deconfining-preconfining phase transition is of second-order kind. 
Large-holonomy calorons dissociate into constituent BPS monopoles, whose masses {Mm) — 
approach zero due to the logarithmic divergence of the electric gauge coupling. Consequently, a 
new ground state emerges as a magnetic monopole condensate. This ground state is responsible 
for a dynamical breaking of the dual U{1) gauge symmetry. The monopole condensate is described 
by an inert complex scalar field (/)j\,f- The gauge mode 7 becomes Meissner massive: 

m^^g{tMM\. (2.105) 

The preconfining character of the magnetic phase is reflected in the fact, that the phase does 
allow for the propagation of massive dual gauge modes, while fundamental, heavy fermions are 
confined by the monopole condensate. 

Due to the genesis of the Meissner mass of 7, a longitudinal polarization arises additional to 
the two transversal polarizations. This causes a discontinuity pc.E PcM in the energy density 
for T \ Tc.E, since the excitation of this additional agree requires energy. 

In the preconfining phase the ground state equation is given as 

p9^ ^ _pgs ^ ttKIjT . (2.106) 

In contrast to the deconfining phase, the negativity of the ground state pressure is not a result of 
monopole-antimonopole attraction, but arises from collapsing and re-created center-vortex loops 
( Abrikosov-Nielsen-Olesen vortices) . 

The confining phase: There is a second phase transition from the magnetic to the center phase. 
The center phase is of confining character. 

If the temperature approaches the critical temperature of the preconfining phase, Tc.m = 
0.83Tc_£;, from above, the magnetic gauge coupling diverges logarithmically at this phase boundary. 
The excitation 7 becomes infinitely massive and a phase transition of Hagedorn type occurs. 

A condensate of massless, single center-vortex loops is generated, due to the decay of the 
monopole condensate. The local magnetic center symmetry is dynamically broken. Massive exci- 
tations are generated by twisting a vortex loop at a self- interaction point. A fascinating property 
of the center phase is that in contrast to the previous phases this vortex loop excitations are spin-i 
fermions. 

In the confining phase the ground state pressure and the ground state energy density are 
evaluated to be precisely zero. 
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Chapter 3 



Correlation of energy density in 
deconfining SU(2) Yang-Mills 
thermodynamics 



The main aim of this chapter is the calculation of the two-point correlation function of the canonical 
energy density (6oo (2^)600 (?/)) a-t different locations x and y and at equal time. This quantity is 
associated with the energy transfer between the two locations x and y. At first, we calculate this 
object in a pure, thermalized U{1) gauge theory, where the propagating gauge field is equivalent 
to the photon of the SM. Our results will serve as a benchmark for the more involved calculation 
of the case when this object is embedded into a deconfining SU(2) Yang-Mills theory, where the 
photon is identified with the massless mode surviving the dynamical gauge symmetry breaking 
SU{2) U{1) in the deconfining phase of SU{2) Yang-Mills thermodynamics. We investigate 
how the energy transfer is affected as mediated by the massless mode interacting with the non- 
trivial massive excitations at low temperatures. Finally, we will propose a possible explanation 
concerning the stability of old, innergalactic, cold clouds of atomic hydrogen. The results presented 
in this chapter are summarized in [39j . 



3.1 General strategy in Coulomb gauge 

The two-point correlation of the energy density is computed by letting derivative operators, asso- 
ciated with the structure of the energy- momentum tensor 6^,^, act on the real-time propagator of 
the C/(l) gauge field. 

Recall that the traceless and symmetric (Belinfante) energy- momentum tensor of a pure C/(l) 
gauge theory is given as 

e^. = -F^^F.x + ^g^.F'^^F.A , (3.1) 

where i^^j, = d^A^ — d^A^j,, and denotes the U (1) gauge field. In the deconfining phase of SU{2) 
Yang-Mills thermodynamics the effective theory contains an inert, adjoint Higgs field (p OIM], and 
the field strength of the abelian theory can in Eq. (|3.ip be replaced by the 't Hooft tensor J^^^, 
to define an SU (2) gauge invariant energy-momentum tensor. This tensor is defined as [26j 

^t^f = TTT'/'aG^iya rn^eabc4>a{D ^(j))b{Dt,(f>)c , (3.2) 

m e\(p\-^ 

where C^^ is the SU{2) field strength of topologically trivial, coarse-grained fluctuations, 
denotes the adjoint covariant derivative, and e is the effective gauge coupling. Obviously, the 
quantity defined by the right-hand side of Eq. p.2p is SU{2) gauge invariant. In unitary gauge 
4>a = 5a'i\4>\ the 't Hooft tensor reduces to the abelian tensor F^^ defined on the massless gauge 
field Al. 
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We evaluate the connected correlation function (6oo(2;)6oo(y)) in four-dimensional Minkowskian 
spacetime (500 = !)• As described in the theory part the energy- momentum tensor 6^*^ is a gauge 
invariant quantity of composite fields. We use Wick's theorem in order to decompose the 2-point- 
function of the composite tensor into a sum of products of 2-point-functions involving the elemen- 
tary gauge fields, i.e. of the propagator Df^^,. Remember that in real-time quantum field theory at 
finite temperature the propagator consists of a vacuum and a thermal contribution [IHl SIl BH US] • 
In Coulomb gauge and momentum space is given as 

D^,{p,T) = D-{p,T) + Dl^,ip,T) 

= -Pj.(p)^^+*^-Pj.(p)2^<5(/)nB(/3bo|), (3.3) 

where 

Pooip) ^ P^M = Pfoip) = 0, 

P^ip) ^ <5.,-^. (3.4) 

In this propagator = (1, 0, 0, 0) can be interpreted as velocity of the heat-bath, /? = ^ denotes 
the inverse temperature and ns {x) = jfzt represents the Bose- Einstein distribution function. By 
virtue of Eq. p.ip one obtains 

(eoo(a;)eoo(2/)> = 2{d^o A^{x)dyO A' {y)){d,oA).{x)dyoAr{y)) 

-goo {d^o {x)dy. A^ {y)){d,o A^ (x)dy^ A^ (y)) 
+500 (9,0 A^ (x)dy^ A^ {y)){d,o A^ {x)dy. A, [y]) 
-goo (5., A^ {x)dyO A^ (y)) {d^. A^ {x)dya Ar{y)) 
+.900 A^ {x)dyO A^ (y)) {d^>. A^ {x)dyO Ar (y)) 

+ ^{d,A\x)dy^A^y)){d,.Ax{x)dy.Ariy)) 
- ^ (9,^ A^ {x)dy^ A^ (y)) {d.,. A^ {x)dy. A,{y)) 
-^{d.^A\x)dy^A^{y)){d,.A^{x)dy.Ar{y)) 

+^{d.. i^)dy. A" {y)) (5,A A^ {x)dy. A^iy)) 

+2{d,oA''{x)dyoA%y)){d.,oAo{x)dyoAo{y)) 
-2{d,oA\x)dy.A''{y)) (9,0 Ao {x)dy^Ao{y)) 
+2{d,rA\x)dy.A\y)){d,^Ao{x)dy^Ao{y)) 
- 2goo {d^«A'>{x)dyoA\y)){d,oAo{x)dyoAo{y)) 
+4500 (9,0 AO {x)dy. AO {y)) {d,o Ao {x)dy^ Ao (y)) 
-2goo{d..^A'>{x)dy.A''{y)){d,^Ao{x)dy^Ao{y)) 

+ ^^{d.,oA\x)dy.A\y)){d,oAo{x)dyoAo{y)) 

-glo {d.o AO {x)dy. A" (y)) {d^a Ao {x)dy^ Ao (y)) 
2 



+ ^{d,^A°{x)dy.A%y)){d,^Ao{x)dy^Aoiy)). 



2 



The last nine lines in Eq. p.5|) arise from the term oc u^ity in the propagator, see Eq. (13. 3p 



(3.5) 



20 




Figure 3.1: Feynman diagram for the correlator {Qoo{x)Qoo(y)) in a pure U(l) gauge theory. 
Crosses denote the insertion of the composite operator Qqq. 

3.2 Two-point correlation of energy density in thermal U(l) 
gauge theory 

In evaluating the expression in Eq. (|3.5p the derivative operators are taken out of the expectation, 
and Eq. (j3.3p is inserted. By momentum conservation the expression Eq. p.5p separates into purely 
thermal and purely vacuum contributions, see Fig. 13. II 

3.2.1 The thermal contribution 

In this section we consider the thermal part of the two-point energy density correlator. We evaluate 
the derivative operators in Eq. (|3.5p and obtain 

(eoo(a;)eoo(2/))"' = 

" J -0)iP^^(p)P=^r{k)poPakok%2nf5{p')5{k^)nB{f3\po\^^^^ 

+ / -0^,P^^iP)PUk)p^P''k^kr{2T,)H{p^)5{e)nBm^^^ 

^ / ^iP'''(p)P^r{k)p'^P^k,k^{2n)H{p^)5{k^)nB{P\^^^^^ 

+ / -0^iP^^(p)P>^r{k)p^Pokxk^{2nf5{p')8{e)nBW\po\)^^^ 

I -§^aP'''^p)p>^-^^)p''^p''^-^^^'^^)^^^p^)^^^^^^ 

" 4^ / -0)iP^'iP)P>^r{k)p^P"kxk.('^T:f5{p^)5{e)nB{P^^^^^ 

(3.6) 
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The tensor contractions are given as 



P--(p)P-(fc) = 1 



Act 

3T «T/„\ dT 



(P ■ k)^ 

p2k2 ' 
(p.k)2 



The Dirac (^-Distribution with quadratic argument can be written in the foUowing way: 



Sip') = 5 (: 



2Vp2 



(3.7) 
(3.8) 
(3.9) 

(3.10) 



This decomposition proves to be useful for the calculation of the integrals concerning the zero- 
coordinate. We use the following property of the i5-distribution: 



(3.11) 



dxf{x)S{x -rf)^ firi) . 
The integration yields four summands: 

{5{pa - Vp^) + 6(po + Vp^)] • [^(fco - Vl?) + 5{ko + Vl?)] = 



= '5(po - Vp')<5(Ao - Vk^) + 6{po - Vp')(5(fco + Vk^) 



-5{po 



V)'5(fco 



6{po 



V)(5(fco + Vk^) 



(3.12) 



Performing the (trivial) integration over the zero-coordinates leads us to 



(eoo(x)eoo(y))"' = 



(pp 



3,p|nB(/5|p|)e'P^ 



(2^) 

d?p f d?k ( pk 



|p||k|nB(/3|p|)nB(/3|k|)e'P^e^''^ 



(3.13) 



where z = x — y. In deriving Eq. (|3.13p we have set = thus neglecting oscillatory terms. 
This prescription should reflect the time-averaged energy transport and is technically much easier 
to handle. An expression for the two-point correlator of energy density including the oscillatory 
terms can be found in Sec. (|B.2p . 

In order to evaluate the integrals, we introduce rescaled momenta and rescaled coordinates 
subject to Pi = (3pi, hi = (3ki and Zi = Zi/ (3, respectively. Subsequently, we introduce 3D spherical 
coordinates 



{Px,Py,Pz) (IpI sin 6'ipi COS 0|p|,|p| sin 6*1^,1 sin (/)|p|,|p| COS 6'ipi) 

{kx,ky,kz) (|k| sinei^cos^,^, |k| sin6'j^sin0^, |k| cos6l,^) , 



(3.14) 
(3.15) 



and express the spatial scalar product between the momenta p and k as function of the polar and 
azimuthal angles 



pk 

|p||k| 



^) + cos 9p cos ( 



(3.16) 



Without constraining on generality, we can set (jjk equal to zero. At first, we evaluate the integration 
over the azimuthal angles (j)p and (jjk- Thus the integrals in the remaining variables |p| and 9 
factorize. 
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We obtain 



(eoo(x)eoo(y))"' = 




(3.17) 



In order to evaluate the ^-integration, we write the spatial Fourier transform in a more con- 
venient way. We are allowed to choose the distance vector z to point into the three direction. 
Thereby, the scalar product of momenta q and distance z in the exponential can be written as 
q • z = |q| |z| cos 9. Finally we use the expansion of a plane wave into spherical harmonics. It holds 
that 

oo 

g^|q||z|cose ^ ^,'(2; + l)j;(|q||z|)P,(cos0) , (3.18) 

where 6 = Z(q, (^), q = p,k, ji denotes a spherical Bessel function of the first kind, and Pi{cos0) 
is a Legendre polynomial in cos 9. An important fact is that the Legendre polynomials obey an 
orthogonality relation defined on the interval (—1, 1): 

+1 

/ dx Prn{x)Pn{x) = —5mn ■ (3.19) 

J 2ni + 1 

-1 

Subsequently, we write the thermal part of the two-point correlator in terms of linear combina- 
tions of Legendre polynomials with the argument cos 9 and exploit the orthogonality relation( |B.9p . 
Consequently, the infinite sum in the plane wave expansion reduces to a finite amount of terms. 
The orthogonality relation of the Legendre polynomials fixes the degree of the involved spherical 
Bessel function. 
We arrive at 



(eoo(x)eoo(y))"' 



(2 



II^D) 



(27r)6^8 



(2 



g /'dipi^^(jo(ipiiz|))-i-^ f d\p\-l^{j2i\p\m 











(3.20) 
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After calculating the binomials and arranging integrals of the same type together, we obtain 



(eoo(x)eoo(y)>' 




d\p\ 



elpl - 1 



Jo(lPl|z| 



rf|p|^#^j2(|p||i 




(3.21) 



We remark, that only partial waves of an even degree contribute to the result, modes with odd 
parity vanish for reasons of symmetry. 

The remaining integral concerning the radial momentum are solvable in an analytical way by 
implementing the following definite integrals [44] : 



dx 



|coth(7r6)-l 



„2m+l 



dx- 



cos{bx) 



- 1 



Q2m+l 



|coth(7r6)-l 



(3.22) 
(3.23) 



where 6 > 0, m > and m e Z is required. Evaluating the last integration, we arrive at the final 
result for the thermal contribution of the two-point correlation of energy density in a U{1) gauge 
theory: 

(eoo(x)eoo(y))"' - 

1 / 647r2 / 1 TT^ coth(7r|z|)cosech^(7r|z| 



(27r)6/38 
327r2 



7r|z|(3coth(7r|z|) + 7r|z|(3 + 27r|z| coth(7r|z|)) cosech^(7r|z|)) 



2 z 



514 



1 



247r4| 



8( - 1 + 7r^|z|^coth(7r|z|)cosech^(7r|i|)) 

+ ( - 8 + 7r|z|(3coth(7r|z|) + 7r|z|(3 + 27r|z| coth(7r|z|) cosech^(7r|z|))))' 



(3.24) 



3.2.2 The vacuum contribution 

In this section we evaluate the quantum contribution to the two-point correlation function 



(eoo(x)eoo(y)) 



(3.25) 



We start again with Eq. p.Sp . but we replace the thermal part of the propagator by the zero- 
temperature quantum propagator. We expect (6oo(2;)0oo(2/))™'' to be a negligible correction to 
(0oo(x)Ooo(y))"' for |z| > 1. and actually the result will show a significant suppression in com- 
parison to the thermal part. Inserting the propagator yields 
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_9oo f d^P f d'^P pXrp P^P" Kkg ,j,._,kc 
2 J (2^)4 J (2^)4 p2 A:2 

7 (2^)4 7 (2^)4^ ' 
_ff00 /" d'^p fj^pXrp P^P" kxkr ,pr_,kC 

2 J (2^)4 7 (2^)4^ ^'''^ p2 p e 

(2^ 7 (^(^0^0 -PoA^oPrA:" +Pa/c Prfc")-^^ ^ 

/■ (i4p /■ d^p ^ e~*P^ e"*'"''' 
7 (2^ 7 J^^Po'^o - '^PokoPrk^ + Pxk Prk'')—^ ^ 

-f /(S^ / Iff*""*" - +f»'>''')i?^ 

(3.26) 



The last three summands are a consequence of the Coulomb- Term in the propagator (|3.3p , which 
describes the 'propagation' of the Aq field. In first instance we use tensor contractions similar 
to Eqs. p.7p . p.Sp . p.9p . Now we accomplish an analytic continuation and rotate to euclidean 
signature (po, ^ Wo, «^o, a;°, 2/° ^ -^2/°, g^iv Sfj.^)- 

We obtain the to Eq. (|3.13p corresponding part resulting from the zero-temperature propagator 



9 f f d'^p 2e'^^\' 1 ( [ d^p . .^e^P^V 

/d'^p f d'^k ( pk \ e*P^ e^'^'' 

^ f d-^p f d^k f pk \ , , ,„ .e*?"; e^'^'^ 

(3.27) 



where C, = x — y and = P/^C/^, fcC = k^(^^,p^ — p^iP^, and k^ — k^k^. The last three lines in in 
the equation above arise from the term oc u^Uv in the propagator, see Eq. (|3.3p . 

In analogy to the calculation concerning the thermal contribution we rescale the variables 
dimcnsionless as p^ = /3p^, and k^ = j3k^. The strategy is now to express the integrals in 
Eq. (I3.27P in terms of four-dimensional hyperspherical coordinates. We restrict to xq = yo to be 
able to compare (8oo(a;)6oo(y))™" with (eoo(x)eoo(y))*''. 
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We define 



P = \p\ 



( cos i/) \ 
sin ■)/) sin d cos 
sin -0 sin sin 
\ sin cos Q j 



< ip,0 <Tr 
with < < 27r 

\p\ = yxf+xf+xfTxf 



(3.28) 



and the integration measure follows as d'^x = |pp sin^ V sin6' dlpldipdOdij)- We proceed now in a 
similar way as in the thermal case. The integration over the azimuthal angles is practicable without 
any difficulties and performed in total analogy to the thermal case. Consequently, the integrals 
over the remaining variables factorize in two identical contributions for each term in Eq. p.27p . 
In order to evaluate the integrals over the first polar angle 9, we determine the spatial distance 
vector z to be collinear to the 3-direction and expand the exponential into partial waves. Due to 
the definition of the 4D spherical coordinates the following expansion is valid 



^i\q\ |z| sin tp cos 



^i'(2? + l)ji(|g||z|sinV')Pi(cos( 



(3.29) 



(=0 



After expressing polynomial factors in cos 9 in terms of a linear combination of Legendre poly- 
nomials Pi(cos6'), we are enabled to exploit the orthogonality relation (|B.9|) . Performing the 
integration over the second polar angle 9 yields 

/ 



(eoo(x)eoo(y))™^ - 



1 



;27r)8/38 



967r^ d\p\\p\'^ dip sin^ ip cos'' 'ipjoi\p\\z\smiJj) 



487r^ 



|p| dip sin ^/icos Vj2(l^ll^l ^^^^) 



967r^ 



d\p\ \p\ dip sin ?/'cosV'ji(|p||z|sin?/') 



327r2 
3 

167r2 



d\p\\p\ dip sin ipjo{\p\\z\sinip) 



IpI'^ dip sin ipj2i\p\\z\sinip) 



d\p\\p\ dip cos ipjQ{\p\\2,\sinip) 



d\p\ \p\ dip sini/'cosi/'Ji(|i5||z| sin?/;) 



247r-' 



d\p\\p\ dip sin ipjo{\p\\z\sinip) 



487r^ 



(^IpI IpI dip sin ipj2{\p\\2\sinip) 



(3.30) 

In Eq. (I3.30p the last four terms result from the term cx u^u^ in the propagator, see Eq. ()3.3p . 
(Details of this straight-forward calculation can be regarded in the appendix). Subsequently, we 
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consider the integration over the second polar angle 6. Merely the terms containing spherical Bessel 
functions of even degree contribute, obviously the third and the seventh term disappear for reason 
of symmetry. 

The integrals can be evaluated by applying Sonine's first integral formula 



7r/2 

j dV' J^,(asinV')sin'^+ii/'cos2''+i V = 2P+^T{p + l)a-P-^ Jp+^+i{a) 



(3.31) 



where denotes an ordinary Bessel function of the first kind, and the generahzation 



r 

diPfiinP^cos''i}Jp{asm^) = 2-^-p{1 + (-l)")a^r ^ ^ 2 ^ ^ ) 



• 1F2 



l + ^ + p 2 + p + u + p 1 + p -a^ 



2 ' 2 ' 2 ' 4 

provided that 9^ei^ > -1, 9lep > -1, lKe(^ + p) > -1 anda > 0). We get 
(Goo (x) 600 (y))™^ 



1 



(27r)8/38 



/ 


/ 00 


967r2 • 


/« 


V 


\o 



+48^2. ( ld\p\--{-QM\p\\z\)~2\p\\z\J2{\m\) 



(3.32) 



4..nn p f'l 3 , (IpII^I)^ 

+3|p||z|iF2 ( 2^2' ' 4 



+ 



+ - 



32^ 
~3~ 

16^2 



d\p\^{3U\p\\i\)-\p\\i\M\pM 



d\p\^M\p\\i\) 



\0 

00 



+ 2ATT-'- I ld\p\^M\p\\i\) 



ASn' ■ - d\p\ 3Jo(|p||z|) + |p||z| Ji(|p||z|) - ^ 1F2 -; 1, -, 



3 A., 3 (|p||z|)^ 



\p\\i\' ^V2' '2' 4 



(3.33) 



The integration over the radial momentum is still remaining. Evaluation of these integrals^ yields 
the final result for the vacuum contribution of the two-point correlation of the energy density in a 



^For a detailed calculation consult Sec. llB.3 
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Figure 3.2: The ratio i?th- 
gauge theory. 



;;u(l) 



(eoo(x)eoo(y))*'' 

(eoo(x)eoo(y)>"'"^ 



as a function of |z| for a therniahzed pure U(lJ 



thermahzed U{1) gauge theory: 



(eoo(x)eoo(y)>^ 




27ry 






F 487r2 • 


167r2 


(-) 


3 





7r4/38 |z 




(3.34) 



The contribution resulting from Coulomb part of the propagator p.3p in Eq. (j3.33p vanishes; 
this observation is in accordance with the fact that no energy transfer between points x and y 
is mediated by the Coulomb part of the photon propagator. In Fig. 13. 21 the ratio i?th-vac;u(i) = 

(6oo(x)eoo(y)) shown as a function of Izl. For example, at T = 5.5, 8.2, 10.9 K a distance 
|z| of 1cm corresponds to |z| ~ 24,36,48, respectively. As a consequence, the thermal part of 
(0oo(x)8oo(y)) dominates the vacuum part by at least a factor of a hundred. 



3.3 Two-point correlation of energy density in deconfining 
thermal SU(2) Yang-Mills thermodynamics 

In this section we consider now the energy transfer mediated by the massless mode surviving the 
dynamical gauge symmetry breaking SU{2) U{1) in the deconfining phase of SU{2) Yang-Mills 
thermodynamics. As already described in sec. ()2.3|) . this symmetry breaking is a consequence of the 
nontrivial thermal ground state composed of interacting calorons and anticalorons. Upon a unique 
spatial coarse-graining this macroscopic ground state is described by a spatially homogeneous field 
configuration, in particular a topological trivial pure-gauge solution a^'^ with an inert, adjoint 
scalar field as static background J^. Due to the adjoint Higgs mechanism Two out of three 
directions in the SU{2) algebra dynamically acquire a temperature dependent mass. Working in 
unitary-Coulomb gauge, where the adjoint Higgs field 0° is given as 0" — 6°'^\4)\ (a — 1,2,3), the 
tree-level massless, coarse-grained, topologically trivial gauge field is A^. 

Since the gauge group is of non-abelian kind, the ground state is not of trivial nature (as in the case 
of an Abelian gauge group like U{1)) and interactions between the massless mode 7 and the two 
massive excitations have to be considered. We investigate how the energy transfer between 
points X and y as mediated by this mode, is affected by this effects induced by the non-trivial 
vacuum . This energy transfer is characterized by the two-point correlator (6oo(x)6oo(y))"' where 
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Figure 3.3: Vanishing irreducible three-loop diagram for the correlator (6oo(2;)0oo(y)) in a ther- 
malized, deconfining SU (2) gauge theory. A wavy (solid) line is associated with the propagator of 
the massless (massive) mode. Crosses denote the insertion of the composite operator 9oo- 

Ooo is now calculated as in Eq. (|3.ip replacing by A^^. We expect that the energy transfer 
in regions of interest for the evolution physics of cold, innergalactic clouds is sizeable suppressed 
in comparison to the counterpart in the standard thermalized C/(l)-gauge theory of quantum 
electrodynamics. 

3.3.1 The thermal contribution 

As already explained in Sec. (|A.4p . the dispersion law of the (tree- level) massless gauge field is 
modified. Since we consider now photon propagation through the non-trivial vacuum of a non- 
abelian theory, we have to take interaction with this medium into account. This non-abclian effects 
manifests in the polarization tensor II^i,. This object is calculated as a summation of one-loop 
self-energies for the on-shell massless mode |451 |46] . see Sec. lA.4l In the effective theory it is 
totally sufficient to account for radiative corrections in the correlator (6oo(x)0oo(y))"' in terms 
of a resummation of the one-loop polarization tensor 11^^ for the massless mode only, see Fig. 13.41 
since higher irreducible loop diagrams yield merely insignificant contributions to the total result 
[17]. According to [46 this polarization tensor leads to a modification of the dispersion law in 
comparison to standard model quantum electrodynamics 

Po = p'^J5o = p' + G(r,|p|,A), (3.35) 

where A denotes the Yang-Mills scale related to the critical temperature Tc for the deconfining- 
preconfining phase transition as Tc = A = A. Photon propagation is then governed by 
Eg. (|3.35p . As described in Sec. (|A.4p . the function G acquires relevance for temperatures not 
much larger than Tc- There is a regime of antiscreening (G < 0) for spatial momenta larger than 
|Paa| 0.2 T. This effect, however, dies off exponentially fast with increasing momenta. For 
momenta smaller than IPhighl ~ 0.1 T and larger than |p,o„| ^ 0.02 T the function G is so strongly 
positive that the propagation of the associated modes is forbidden (total screening). 

Initially we derive the photon propagator in deconfining SU (2) Yang- Mills thermodynamics. 
Considering the modification in the dispersion law, Eq. (|3.35p , the propagator follows as 

D,^(p, T) = -Pl(j>) ( p2_G(r,|p|,A)+»6 + '^''^^P' IpI' A))"s(/3bo|)) . 

^This can bo made manifestly SU{2) gauge invariant by substituting the 't Hooft tensor '26' for the field strength 
into Eq. f3A) . 
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Figure 3.4: Feynman diagram for the correlator (6oo(2;)6oo(?/)) in a thermalized, deconfining 
SU{2) gauge theory including radiative corrections to lowest order. (Blobs signal a resummation 
of the one- loop polarization for the massless mode.) Crosses denote the insertion of the composite 
operator 800 ■ 

(3.36) 

After inserting the thermal part of Eq. p.36p into Eq. p.Sp and implanting the same tensor contrac- 
tions Eqs. p. 7113. 8113. 9p like in the U{1) case, we decompose the (5-distribution similar to Eq. p.l2p . 
Performing the integration over the zero-coordinate, the thermal part (8oo(x)8oo(y))"' calculates 
in analogy to Eq. (|3.13p as 




«i3(/3Vp' + G)ns(/3v/k2 + G) c'P(--y)c''^(''-y) 

+2 7(2^7 (2^ (vWJ TF+gTP+I'' 

ns(/?v/p' + G)ns(/3v/k2 + G) e'p(— y)e^k(x-y) 

(3.37) 

The integration over the p^- and /cQ-coordinates is executed without constraints and just fixes the 
dispersion law p.35p . We proceed now in exactly the same manner as in the U{\) case, since 
we introduce spherical coordinates and evaluate the integrals over the angles analytically. Upon 
performing the azimuthal integrations for each summand in Eq. (|3.37p the respective expression 
reduces to a square of two integrals over the polar angle and the modulus of spatial momentum, 
respectively. The integral concerning the polar angle is computed by expanding the exponential 
in spherical harmonics, writing polynomial factors in terms of linear combinations of Legendre 
polynomials and exploiting their orthogonality relation, Eq. (|B.9|) . The calculation of integral 
over the modulus of the spatial momentum is a more complex task, since we have to consider 
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the screening effects induced by the one-loop polarization tensor. We continue the evaluation in 
analogy to the derivation of the black-body spectrum in [IS]- Employing the modified dispersion 
law, the integral over momentum- modulus is replaced in favor of an integral over frequency w. 
Those values of uj, which yield an imaginary modulus of the spatial momentum due to strong 
screening, are excluded from the domain of integration. This regime of strong screening is in the 
range lui < lo < uj2, where uji , 0J2 denote the solutions of the equation 



uj^ - G(tj,r, A) = 0. 



(3.38) 



Finally, we introduce a dimensionless frequency Q = (3u) and a dimensionless screening function 
G = 0^G and arrive at 
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(3.39) 



The frequency integral is calculated numerically. The results for different temperatures and con- 
clusive implications are analyzed in Sec. (|3.4p . 



3.3.2 Estimate for the vacuum contribution 

A this point we would like to obtain an order-of-magnitude estimate for the vacuum part of the two- 
point correlation of the energy density concerning the massless mode in deconfining SU(2) Yang- 
Mills thermodynamics. To do this we ignore the modification of the dispersion law in Eq. p.35p . 
Anyways, the function G has so far only been computed for external momentum p with — 0. The 
difference as compared to the U{1) case is then a restriction of the (euclidean) four- momentum 
p as p2 < 02 existence of a scale of maximal resolution |0| in the effective theory. 

This compositeness constraint , Eq. (j2.96p originates from spatial coarse-graining. Again, we will 
see that (6oo(x)6oo(y))™'' is a negligible correction to (0oo(x)6oo(y))"' for physically interesting 
distances. 

In case of (6oo(x)0oo(y))™° we restrict to xo = yo to be able to compare with (0oo(x)6oo(y))"'- 
Introducing the dimensionless modulus of as = /3(f>, proceeding in a way analogous to the 
derivation of Eq. p.30p . and performing the ip- and |p|-integration^, we arrive at 



^Concerning a detailed calculation of the integrals we refer to the appendix Sec. lB.3l 
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(3.40) 



where 1-F2 is a hypergeometric function and Jq, Ji, J2 are Bessel functions of the first kind (con- 
ventions as in [H]). 



3.4 Numerical results and discussion 

According to sec. p.2p it is evident that if photon propagation is governed by a U{1) gauge theory, 
neglecting the vacuum contribution in comparison to the thermal part of the two-point correlation 
of energy density is well justified. In the following we turn to the case of energy transfer mediated 
by the massless gauge boson in deconfining Yang-Mills thermodynamics. 

The frequency integral in Eq. (|3.39p is evaluated numerically. We compare the estimate for 
the vacuum contribution in Eq. (|3.40p with the thermal part of the correlator (0oo(x)8oo(y))- To 
make contact with SU {2)cmb, whose Yang-Mills scale is A = 1.065 x 10~*eV (13 US], we relate at 
a given temperature the dimensionless distance |z| to the physical distance in centimeter^. 

We define 

(9oo(x)9oo(y))-'- 

''™^'^'^=l(eoo(x)eoo(y)>-r ^'-''^ 

where we use the estimate in Eq. p.42p for |(6oo(x)6oo(y))™°|- In Fig. [3?5lthe quantity i?th-vac;su(2)(|z|) 
is depicted for various temperatures specializing to the case of SU (2) cmb- Notice the strong dom- 
inance of the thermal part. Notice also that although this result resembles qualitatively the result 
of Fig. 13. 21 for fixed distance and varying temperature this is not true for fixed temperature and 
varying distance. Namely, the existence of a nontrivial thermal ground state in deconfining SU{2) 
Yang-Mills thermodynamics constrains quantum fluctuations of the massless mode to be softer 
than the compositeness scale Owing to its trivial ground state, no such constraint exists in a 
thermalized U{1) gauge theory. As a consequence and in accord with Fig. 13. 21 quantum fluctuations 
dominate thermal fluctuations at small distances in such a theory. 



*We can convert the values from Planck units to SI units by applying the relation 

hc=l = 197Mevfm (3.41) 
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3.0 X 10*' 
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Figure 3.5: The function i?th-vac;su(2) (|z|), defined as in Eq. (|3.42p . when specializing to the case 
of SU(2)cMB {Tc — 2.73 K), for various temperatures: black curve (T = 1.5 Tc), dark grey curve 
(T = 2.0 Tc), grey curve {T = 2.5 TJ, light grey curve (T = 3.0 T^). The dashed line separates 
distances smaller than |(/)|^^ from those that are larger than 101^^. 



By virtue of the results shown in Fig. 13. 21 and Fig. 13. 51 we neglect the vacuum contribution to 
(600 (x) 900 (y)) in the following. 

Let us now turn to the interesting question of how much suppression there is in the two-point 
correlation of the photon energy density in the case of SU (2)cmb as compared to the conventional 



U{1) case. In Fig. 13. 61 the ratio R, 



^•th,SU(2)— th,U(l) 



Izl), defined as 



„ _ (eoo(x)9oo(y))-'-^"'^' 



(3.43) 



is depicted for various temperatures as a function of distance in centimeters. Notice the significant 
suppression of the correlation between photon energy densities in the case of an underlying SU{2) 
gauge symmetry as compared to the case of the conventional C/(l) up to a factor of two. Since 
the two-point correlator (9oo(x)6oo(y)) is a measure for the energy transfer between the spatial 
points X and y in thermal equilibrium and hence a measure for the interaction of microscopic 
objects emitting and absorbing radiation, we conclude that this interaction is, as compared to 
the conventional theory, suppressed on distances ~ 1 cm if photon propagation is subject to an 
SU{2) gauge principle. To make the situation even more explicit we have computed the Coulomb 
potential V{r) (r = |x|) of a heavy point charge in case of photons being described by a U{1) and 
an SU{2) gauge theory. This potential is given in a U{1) theory as 
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sin^ 



(3.44) 



Going from U{1) to SU{2), we take into account the resummed one-loop polarization [46^ by 
letting + G in the denominator of the integrand in Eq. p.44p . Here G is the same 

function as discussed in Sec. (|A.4p . The physical situation is a heavy point charge immersed into 
the SU (2) plasma with the photons associated with it being (anti)screened by nonabelian, thermal 
fluctuations. Although the function G is known for on-shcU photons — p^ only this recipe 
should work well for sufhciently large distances since the off-shcllness of the photon, which does 
not mediate any energy transfer from the source, then is sufficiently small. Thus for SU (2) we 
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Rth;SU(2)-th;U(l) 
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Figure 3.6: The function Ra^ su{2)-th.u{i){\'^\), as defined in Eg. (|3.43p . when speciahzing to the 
case of SU(2)cMB (Tc = 2.73 K), for various temperatures: black curve (T = 1.5 Tc), dark grey 
(T = 2.0 Tc), grey (T = 2.5 T^), light grey (T = 3.0 T^). Notice the regime of antiscreening for 
small \z\. For T = 3.0 Tc, where the approximation ~ deviates sizably from the full result for 
the screening function G, we have used G as obtained 1-loop selfconsistently ^5] . 
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1 
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(3.45) 



where the last integral is performed numerically. In Fig. l3.7l both potentials, Vu(i) ('') and V'su(2)cmb 
as well as their ratio are plotted as functions of r. Notice that the dimensionless quantity 



SU(2)CN 



^{r)/Vv(i){r) shows similar suppression as the dimensionless quantity i?th,su 



(2)-th,U(l)(|z|) 

depicted in Fig. 13. 61 This seems to confirm the validity of our above approximation. 

The following tabular lists the maximal resolution \(f>\~^ as dimensionless quantity as well as 
the conversion into the physical interesting unit [cm] for various temperatures. 



Dimensionless temperature A 


\4>\-' 


|(/)| [cm] 


1.5xAc==20,81 


15.10 


0.84 


2.0xAc==27.74 


23.25 


0.98 


2.5xAc=34.68 


32.50 


1.09 


3.0xAc=41.61 


42.72 


1.19 


3.5xAc=48.55 


53.83 


1.28 



3.5 Stability of clouds of atomic hydrogen in the Milky Way 

A large fraction of the interstellar medium in our galaxy consists of warm 10^ — W^K) and 
cold (~ 50 — lOOK) atomic hydrogen gas. The standard theory of the interstellar medium implies 
that the colder the region, the higher the affinity of atomic hydrogen to form molecular hydrogen 
This leads to the assumption that extrem cold regions are dominated by dense clouds of molecular 
hydrogen. 

Initially, we should give a short review concerning the evolution of such cloud structures [8 . In 
the galactic center hot ionized hydrogen gas is emitted. This gas escapes the main galactic disk and 



®A point charge immersed into the plasma locally distorts the latter's ground state, and, strictly speaking, the 
theory to describe this distortion yet needs to be worked out. But measuring the (anti)screening of the potential 
sufficiently far away from the location of the charge should still be describable in terms of unadulterated SU(2) 
Yang-Mills thermodynamics. 
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Figure 3.7: Left panel: Plot of the potentials Vu(i)(r) (dashed line) and 14u(2)cmb (^) (solid line) 
as a function of r at T = 2.5rcMB ~ 6.8 K. Right panel: Plot of the ratio Vsu(2)cmb(^)/^u(i)(?') as 
a function of r. The temperature is set to T = 1.5 (black), T = 2.0 (dark grey), T = 2.5 
(grey), and T = 3.0 Tc (light grey). For T — 3.0 Tc, where the approximation = deviates sizably 
from the full result for the screening function G, we have used G as obtained 1-loop selfconsistently 
PU] , Notice how close to unity the curve at T = 3.0 Tc is as compared to the curve at T — 2.5 Tc 
for the right-hand side panel. 

forms a giant 'mushroom' structure above it. In the course of time the ion gas expands, governed 
by gravitational effects, cools down, the ions recombine and trickle down onto the galactic disk in 
to the voids between the spiral arms of the galaxy, where clouds of atomic hydrogens are formed. 
The clouds continue to cool down to 10 — 25K. At the same time neutral atoms are assumed to 
combine to molecular hydrogen HI via dipole interaction (i.e. SM QED). Numerical simulations 
of the cloud evolution predict a time scale of less than lO^y [5D1 IM] for the process of generating 
a substantial fraction of molecular hydrogen. 

But in 2001, the discovery of an arc-like structure located in-between the spiral arms of the 
Milky Way and consisting of atomic hydrogen HI, designated as GSH139-03-69, was reported in 
[IlIH]. This object GSH139-03-69 is impressive, with a size ~ 2kpc, a mass of about 1.9 • 10^ solar 
masses, and an atomic number density of ~ 1.5cto~^. The cold cloud exhibits a mean brightness 
temperature Tb 20A' with cold regions of Tb 5 .. . lOK. Although GSH139-03-69 is as cold 
as a molecular cloud, it consists mostly of atomic hydrogen. Regarding the observed data, one 
expects that the age of GSH139-03-69 is at most lO^y. However, the astonishing fact is that this 
object possesses an estimated age of about 5 • lO^y, which is an obvious contradiction to the SM 
of the interstellar medium and hence not conform with the SM of particle physics. 

The results obtained in Sec. p.3p concerning the electromagnetic interaction of microscopic 
objects yield immediate implications for the gradual metamorphosis of cold (T = 5 . . . 10 K) astro- 
physical objects such as the hydrogen cloud GSH139-03-69. Although a quantitative estimate of 
the increased stability of atomic hydrogen clouds due to the SU{2) effects in thermalized photon 
propagation is beyond the scope of the present work, Fig. l3.6l clearlv expresses that the mean energy 
transfer in the photon gas of T = 5 • • • 10 K is suppressed by up to a factor of two at interatomic 
distances in the hydrogen cloud GSH139-03-69. It would be interesting to see how the simulation 
of the cloud evolution, taking into account the effects as expressed in Fig. 13. 61 would increase the 
estimate for its age as compared to the standard picture. 

As already pointed out in [4^, the propagation of the 21-cm line, which thermalizes structures 
such as GSH139-03-69, is unscreened even when subjecting photons to S'C/(2)cmb- A recent full 
calculation of G shows, however, that this result is an artefact of the approximation = 0: 
Thermalization of the hydrogen cloud thus takes place solely via the coupling of the photon to the 
nontrivial thermal ground state [52] . 
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Chapter 4 



The spatial string tension in 
deconfining SU(2) Yang Mills 
Thermodynamics 



4.1 The Wegner- Wilson loop and magnetic screening effects 

In this section we investigate a gauge-invariant construction called Wegner- Wilson loop. As al- 
ready mentioned in the Introduction the gauge theories in their strongly-coupled regime are not 
considerable in the framework of perturbation theory, standard Feynman diagrammatic techniques 
are not applicable. Thus in the second half of the 70's the investigation and description of gauge 
theories with a large coupling (for instance QCD in the infrared) demanded the development of 
nonperturbative methods involving for example the Wegner- Wilson loop. This quantity is defined 
as the expectation value of a path-ordered exponential of a gauge field 



where the contour integral is evaluated along a rectangular path around the origin. Gauge- 
invariance is guaranteed by the cyclicity of the trace. 

This object was originally motivated in lattice gauge theories as an attempt at a nonperturbative 
formulation of QCD in its strongly-coupled region. Nevertheless, the asymptotic behavior of the 
Wegner- Wilson loop is an appropriate criterion to distinguish between the phases of deconfinenient 
and confinement. We focus on another aspect of the Wegner- Wilson loop which monitores screening 
effects in a plasma originating from the existence of magnetic quasi-particles (monopoles) at finite 
temperature. 

Consider a hot plasma in the framework of an abelian gauge theory. It is a well-known fact, that 
the existence of long-range static magnetic fields is possible, and no magnetic screening effects oc- 
cur [20]. In contrast to these long-range electric fields are absent due to Debye screening. However, 
in the non-abelian scenario of an SU (N) gauge theory, apart from Debye screening in the electric 
sector, additional screening effects in the magnetic sector emerge. This observation suggests the 
existence of magnetic monopoles in the plasma. This is mathematically corroborated in the spa- 
tial Wegner- Wilson loop, which measures the magnetic flux induced by magnetic quasi-particles 
through the enclosed area of the loop, due to Stokes' theorem. A result from lattice simulations is 
that the thermal average of the Wegner- Wilson loop obeys an area-law, which corresponds to a gas 
of free screened magnetic monopoles. For more details concerning the screening effects resulting 
from magnetic activity due to monopoles at high temperature we refer to |531 154] . 

The spatial string tension is defined as the negative logarithm of the thermal average of the 
Wegner- Wilson loop divided by the enclosed area A of the loop: 




(4.1) 



ln(tr7' exp ie dz^A^) 



(4.2) 



A 
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Figure 4.1: Illustration of the diagrammatic approach. Every summand in this formula represents 
a whole class of A^-gluons exchange diagrams, endowed with an combinatoric factor. 



where denotes the gauge field and e the (effective) coupling constant. The integration contour 
C is given as the edge of a square with sidelength L in the limit L — > oo. 

We are interested in the evaluation of the spatial string tension in the deconfining phase of 
SU{2) Yang-Mills thermodynamics using the effective theory [9]. In this context, we assume 
magnetic monopoles liberated from dissociated calorons to be responsible for an area-law. 

If the spatial string tension shows area-law behavior, we will obtain the following relation (for 
r>Tc): 



^lim ^IntrPexp ^ie j> dz^" = F{T) cx 



(4.3) 



Otherwise, if the spatial string tension is governed by an perimeter-law, we will get 

lim - Intr-Pexp I ie i, dz^^A^" I = /(T) cx (4.4) 



; jdzf'A'^ = f{T) cx 



L — >oo L 

In the following the Wegner- Wilson loop will be denoted as W[C] 

4.2 Calculation of the spatial string tension including the 
on-shell polarization effects 

In order to calculate the spatial string tension, we use an expansion into loops in the effective theory 
|55j . A resummation of all possibly contributing A'^-gluon exchange diagrams can be rewritten in 
form of an exponential of the one-gluon exchange diagram. This fact is illustrated in Fig. 14. II Thus 
we can write for the logarithm of W[C] 

\iiW[C] = -]^Cf j>dx^dy,D^,{x-y), (4.5) 

where Cp denotes the Dynkin index in the fundamental representation, defined as the normalization 
factor of the generators of the algebra (i.e. Cp = \)^ and 

3 

^M^-E4"^ (4.6) 

a=l 

is the sum of the tree-level propagators (|A.36[rS~34|) of the three gauge modes in deconfining, effec- 
tive Yang-Mills thermodynamics. On-shell polarization effects concerning the tree-level massless 
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gauge mode are taken into account by including the resummed one-loop polarization tensor into 
the propagator of the massless gauge mode (IA.34P . We constrain ourselves to the magnetic part of 
the polarization tensor (screening function G(p°,p)), electric sreening effects (screening function 
G{p°,p)) do not concern the spatial Wegner- Wilson loop. We have to emphasize that we only 
consider on-shell polarization effects for the massless gauge mode, radiative corrections for the 
massive gauge modes are not taken into account. Beyond this fact we should remark that this 
on-shell polarization tensor for the massless mode is calculated exactly up to one- loop order, since 
higher loop contributions vanish identically in the effec five theory [37], see Sec. lA.41 

The tree-level propagators in position space are given as the Fourier transformations of their 
momentum space counterparts 
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where the transversal projection operator is given as 
Poo(p) = Pl{p)^Pl{p) = 



(4.8) 



(4.9) 



We calculate the contour integral in the 1-2-plane, that is a;o = yo = a^o = J/3 = 0, and consider 
at first an arbitrary tensor structure D^^{p). Later we will insert the explicit tensor structure of 
D^'J for the massive gauge modes {V^ gauge modes) and 1?^^^ for the massless one (7 mode). We 
obtain 
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We arrive at 
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At this point we insert the propagators for the and th 7 mode. In case of the massive gauge 
modes we only have to consider the thermal part of Eg. I|A.36[) . since due to the constraint on 
the maximal off shellness Eq. ()2.96p the vacuum part is vanishing. Subsequently, we perform the 
integration over the po-coordinate in the thermal contributions. Thus we obtain 
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We rescale the momenta and the squared gauge boson mass to dimensionless variables in terms 
of temperature and the loop length: 



Pi ^Pi- L 
and with e = \/8n 

„ 2 m'^ _ (2e)2 A 
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T2 T2 27rT A3 " 
The screening function G (as dynamically generated mass) is rescaled in the same manner 

G- ^ 



(4.12) 



(4.13) 



(4.14) 



Subsequently, in order to control the limit of growing loop length L ^ 00, we introduce the 
dimensionless parameter t into the integral (|4.1ip via the definition 



T = T-L. 



(4.15) 
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Figure 4.2: Two-loop contribution to the pressure. Cutting of the massless hue yields exactly 
on-shell one-loop polarization tensor for the massless gauge mode. 



Finally, Eq. ()4.1ip is recast as 
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(4.16) 

An area-law behavior of the spatial string tension corresponds to the existence of a nonvanishing 
limit 



= — lim 



lnT4^[q 



(4.17) 



a perimeter-law is given, if a finite limit according to 
..-hmi^.T^ 



(4.18) 



exists. Let us analyze the asymptotic behavior of Eq. (|4.16p at high temperature and low momenta 
for T — > 00: It is obvious, that in this limit the only contribution results from the massless mode, 
the contribution of the massive modes vanishes. This can be seen by an expansion of the integrand 
of In ^^[CJv'i, in particular the Bose-Einstcin distribution, for the momenta. In the limit t 00 
the integral is proportional r^^. Concerning the thermal contribution of the massless mode we 
observe a perimeter-law for each temperature in the limit r — > 00. However, for values of t below 
T*, determined such that L ■ we are able to observe curvature in lnW^[C]I^''. Concerning the 

vacuum part the integrand can be estimated by neglecting the screening effects. We observe that 
lnl/F[C]™'' represents a r-independent, infinite constant, and thus does not contribute to the force 
associated with the potential of a fundamental quark and its antiquark in the (2-|-l)-dimensional 
situation described by the spatial Wilson loop. A detailed discussion is beyond the scope of this 
thesis and will be part of future work 56J. 

The negative logarithm of the Wegner- Wilson-loop W[C] is depicted as function of the param- 
eter T in Fig. 14. 31 Including merely the on-shcU polarization tensor for the massless mode, we are 
obviously not able to reproduce the generally accepted results from lattice simulations. This leads 
to the idea that the lattice-obtained high-temperature behavior of the string tension is a conse- 
quence of off-shell contributions to the polarization tensor. As mentioned above, an area-law of the 
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-In W[C]y 




Figure 4.3: Plots of — lnT4^[C]!y'' as a function of r = T-L. The dotted lines correspond to the value 
of L coinciding with the minimal length scale \4>\^^ in the effective theory at a given temperature. 
The left hne is for T = 2Tc, the right line for T = STc, and the hnc for T = AT^ would be at 
T*[4rc] = 65.77 and thus is not contained in the figure. 

string tension requires the magnetic activity of free screened magnetic monopoles. These isolated 
magnetic charges are liberated during the dissociation of calorons. The likelihood for the decay of a 
large-holonomy caloron into a monopole-antimonopole pair is enormously Boltzmann-suppressed, 
^g-/3M„/„ ^ c^^'^ ^ . However, isolated magnetic monopoles and antimonopoles can be generated 
from small-holonomy calorons by dissociation. A small-holonomy caloron interacts with an off-shell 
TLM-mode. In this process an amount of energy Ai? is transferred to the bound system of th e 
monopole-antimonopole pair. If llE is large enough to exceed the attractive potential induced by 
quantum fluctuations [S^, the pair of constituents can split up into isolated magnetic charges. 

Another aspect of this scenario concerns the two-loop contribution to the pressure AP^^ . This 
quantity was computed in [3^1 ET] . After scattering at the caloron, the TLM-mode yielding the 
energy to decouple the monopole-antimonopole pair possesses diminished energy than in the initial 
situation. As a consequence, the thermodynamical pressure is diminished and isolated magnetic 
monopoles are liberated. Since the pressure is proportional to — T'*, we can conclude that the 
density of isolated magnetic monopoles grows with increasing temperature. Consider the two-loop 
Feynman diagram concerning the pressure, where the massless line is cut. As depicted in Fig. 14. 21 
we obtain exactly the polarization tensor for the massless mode. Thus we can interpret the effects 
summarized in the one-loop polarization tensor on a microscopic level as scattering of massless 
modes on calorons. This leads to isolated magneti c charges and thus to magnetic flux in the 
plasma, measured by the spatial string tension. Finally, we have to emphasize that the off-shell 
polarization tensor contains an additional contribution (resulting from diagram A, see Fig. lA.3|, 
in contrast to its on-shell counterpart. This fact let us assume that in our effective theory a 
reproduction of the area-law behavior of the spatial string tension at high temperature possibly 
requires a consideration of the off-shell polarization tensor. 
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Chapter 5 

Summary and Outlook 



In this thesis wc have considered two aspects of SU{2) Yang-Mills thermodynamics in the decon- 
fining phase. At first, we have calculated the two-point correlation function for the energy density 
of the photon, followed by an evaluation of the spatial string tension including on-shell polariza- 
tion effects for the massless mode in the effective theory. This investigation was motivated by the 
postulate that the U gauge group of the SM is not the progenitor to the gauge group of QED. 
Alternatively, the U{1)y gauge group is originating from a more fundamental SU{2) gauge group, 
named S'J7(2)cmb, which dynamically breaks down to C/(l)y. 

In this effective thciory the thermal ground state is composed of interacting calorons and an- 
ticalorons with trivial holonomy and characterized by an adjoint scalar field </> and a pure-gauge 
configuration aff . The macroscopic scalar field </> emerges as a result of spatial coarse-graining and 
acts as a Higgs field, responsible for the dynamical symmetry breaking, with two out of the three 
gauge modes acquiring a temperature dependent mass. The remaining massless gauge boson is 
identified with the photon. Distorted photon propagation according to the postulated SU{2) gauge 
principle manifests itself in the polarization effects due to interactions with the massive excitations. 
As a consequence, the dispersion law is modified. 

In Chapter three we have calculated the two-point correlator of a composite object, the energy 
density of the photon, in the conventional U{1) gauge theory of electromagnetism and in deconfining 
SU (2) Yang-Mills thermodynamics. This correlation function can be interpreted as a measure for 
the energy transfer between two spatial points. Apart from the fact that the computation is 
interesting from a technical point of view, we have considered a new feature of photon propagation 
governed by the SU{2)cmb gauge theory. The visible suppression of (6oo(x)6oo(y)) in case of 
photon propagation underlying a SU{2) gauge theory in comparison to the conventional SM-QED 
implies a possible explanation for the unexpected stability of large, cold, innergalactic clouds 
consisting of atomic hydrogen. 

Subsequently, we have considered the spatial string tension in deconfining SU{2) Yang- Mills 
thermodynamics. This quantity can be regarded as a measure for the magnetic flux induced 
by magnetic monopoles at finite temperature through the area enclosed by a spatial contour. 
Including only on-shell polarization effects for the photon, we were not able to reproduce the 
lattice-obtained area-law. This has led us to the conclusion that the area-law is a consequence of 
off-shell contributions to the one-loop polarization tensor, which is natural since the Coulomb-field 
of a monopole is composed of oflf-shell modes. A calculation of the off-shell radiative corrections 
and thus a possible verification of our statement about the origin of the generally accepted area-law 
for the spatial string tension would be an interesting challenge for future works. 
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Appendix A 

Thermal field theory 



In this section we summarize the fundamental concepts of finite temperature field theory. The 
analogy of quantum statistical mechanics and quantum field theory is developed, in order to transfer 
methods of zero-temperature quantum field theory to thermal field theory. The imaginary-time 
and the real-time formalism in thermal field theory are reviewed, and the free propagator in a 
thermalized theory is derived. After that, the Feynman rules in QFT at finite temperature are 
stated, which are necessary prerequisites for the evaluation of Feynman diagrams. For a detailed 
introduction into the topic of finite temperature field theory and its applications in high energy 
physics, we refer to the literature [iO l lil ] li^ l . 



A.l The partition function as functional integral 

We restrict ourselves to the bosonic case. 

At first, we remind the connection of quantum statistical mechanics and quantum field theory. 
Consider the amplitude A{(j)f ,tt] (j)i,ti) for the evolution of a system from an initial Heisenberg 
state at time U to a final state \(j)f) at time tf. The dynamics of the system is governed by 
the time-independent Hamiltonian 

H = Jd^xni(l),n) ^ J d^x{TTdo(t)- C{(j),7r)) , (A.l) 

where Ti. and C denote the Hamiltonian density, respectively the Lagrangian density, as a func- 
tional of the field cf) and its conjugated momentum tt. If we introduce the time-evolution operator 
Q-iH{t}-u) ^ we are able to write the transition amplitude as 

A{<j,f^uc|>^.U) = (0/|e-^^(*/-*')|0.). (A.2) 

In the path integral formulation of quantum field theory this object is calculated via functional 
integration in phase space: 



A{(t>f,tf;(t>,,U) = f P0 /ptt exp I i /di /d^a; ( 7r(x,t)^^^l^ -H(0(x,t),7r(x,<)) 



(x,,t,) 



(A.3) 



In the following, we assume that Ti, depends only quadratic on the conjugated momentum. In this 
case one can complete the square in the exponential, thus the 7r-integral is of gaussian type and 
can be evaluated analytically. This leads us to 

A{cj)f,tf-cj),,U)^Af J l?0e^^[*l, (A.4) 
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with the action functional = J d^x C{(j)) and a normalization constant A/" (which we can absorb 
in the integration measure). The integration considers every possible field configuration on which 
the system can evolve from the initial to the final state, weighed with a phase factor proportional 
to the action of the field configuration. 

In several cases it is convenient to evaluate path integrals in euclidean instead of Minkowskian 
spacetime. Especially in quantum field theory at non-zero temperature euclidean path integration 
is a common method of calculation. We obtain the path integral in euclidean formulation by 
performing an analytical continuation from real to imaginary time, 

t —ir, or —1x4 , (A. 5) 

where r real. This transformation is called Wick-rotation and changes the signature of our metric. 
Since in Minkowskian signature we have x^Xf^ = — we obtain x^x^j, = x^x^^ = -f- in 
euclidean signature. The transition amplitude for the field is then (/)(x, t) ~ (/)(x, —ir) is then 
transformed into 

A(^/,T/;</)„T,) = (.^/|e-^(^/-^')|)=AA [Vcpe-^^^^K (A.6) 



Here S'£;[(/)] denotes the action functional in euclidean spacetime, the boundary conditions concern- 
ing the paths are given as 4>i{x) — 0(x, r^) and 4'f{'x-) — (/"(x, Tf). (Without constraining generality, 
we set Ti = ti = 0.) 

Now we look at quantum statistics in the canonical ensemble. Our intention is to transfer the 
concepts of functional integration in quantum field theory to statistical mechanics with the aim to 
develop a path integral formulation of the partition function. The partition function occupies a 
central role in the canonical ensemble of statistical mechanics and is given as 

Z = trc-'3^ = ^(0„|e-^'«|0„). (A.7) 

n 

The sum is to evaluate over a complete set of states of the system with Hamiltonian H and inverse 
temperature /? = ^. The similarity of the partition function and the euclidean path integral 
formulation of the transition amplitude is therefore obvious. Formally, we are allowed to write 

Z = ^^((/.„,-^r;0„,O). (A.8) 

n 

If we compare this formula with the euclidean expression for the transition amplitude, we can 
write the canonical partition function as an functional integral, since we adjust the following 
requirements. We identify finite temperature with imaginary time compactified on a circle with 
circumference (3 and interpret e~^^ as imaginary time evolution operator. As consequence of the 
trace operator in the partition function, the physical fields have to satisfy /?-periodic boundary 
conditions in euclidean time, the functional integration is constrained to periodic fields, 0(x, 0) — 
(j){x,i3). With these conditions we obtain 

Z= J P^e-^^I"^!. (A.9) 

</.(x,O) = 0(x,/3) 



A. 2 Derivation of the propagator for the free scalar field 

Now we receive the propagator by applying a standard functional method in quantum field theory. 
We define a generating functional Z{(3, J) with a source term J through 




-SE{P) + !id''xJ(x)<j,(x) _ 



It is obvious that for a vanishing source term, the generating functional is equal to the canonical 
partition function Eq. (IA.9p . 



44 



Functional differentiation of Eq. (|A.10|) yields the propagator in euclidean, that is imaginary 
time: 



1 S^ZiP,J) 



Z{P) 5J{x)5J{y) 



1 



,=0 Z{I3) 



(A.ll) 



We specialize to the case of a free bosonic field. The free generating functional in euclidean 
formulation is the given as 







ZfW.J) 



I?0 exp 



dr / d X (h(x 



1 



(x) — J(x)0(x) 



= ZfW) exp \^ J d^xd^yJ{x)^F{x-y)J{y) 
The two-point Greens function /S.p{x — y) solves the partial differential equation 



(A.12) 



9r2 



V' + m'\ (x =<5(r,-T^)<53(x-y). (A.14) 



As solution in Fourier space, we obtain the imaginary-time propagator, also called Matsubara 
propagator 



1 



(A.15) 



where the aj„ denote the discrete Matsubara frequencies ujn 



and ujk is defined as ujk 



-\- . The corresponding expression in position space is given as 



d^k 



/3 ^ y (27r)3 w2 + /j2 + ^2 ■ 

A comparison with the zero-temperature counterparts in euclidean metric 

A{x-y) = 



d'^k e^''(^-y) 



(A.16) 

(A.17) 
(A.18) 



(27r)4 fc2 + ' 

where k = (fc4, h) , k^ = kf + , yields transformation rules from zero to finite temperature: 

ki UJn, (A.19) 



2tt 



/(fc4 



d^X 



dr / d^x . 



(A.21) 



The temperature dependence is contained in the infinite sums over the temperature dependent, dis- 
crete Matsubara frequencies. As a consequence of the frequency summation, the finite temperature 
propagator A^^\x — y) consists of a zero-temperature and a finite-temperature contribution. 

A frequency sum -g f{uJn,K), where conventionally K represents the remaining variables 
on which / depends on, can be executed with methods of complex analysis. We substitute the 



^We should mention, that 5{tx — Ty) represents a periodic <5-Distribution 



(A.13) 
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Figure A.l: Integration contour. 



discrete values ujn by a continuous variable uj and define a function f(Lu,K), which we assume to 
have no singularities on the real axis. Additionally, we define the complex function 

h{Lo) = ■ = il3 , (A.22) 

whose poles are located at a; = with unit residue. Employing the residue theorem, we are able 
to write the sum as follows, 

F{K) = ^jdLo h{cu)fiuj, K) , (A.23) 

a 

where a designates the contour shown in Fig. lA.ll With respect to the behavior of h{uj) for large 
imaginary parts, Eq. (jA.23p can be decomposed into three integrals 

00 00 — ze oo+^e 

F(K) = ^ [ d^f{u:,K) + ^ [ d^ip^ + ^ [ (A.24) 
^ ' 2t: J ■'^ ' ' 271 J e'P'^ - 1 27r 7 e-^'^'^ - 1 ^ ' 

— 00 — 00— — 00+ie 

If we provide that f{uj,K) is meromorphic as a function of lu and well-behaved for — > 00, 
the integration contours concerning the second and third integrals can be closed in the lower, 
respectively upper half plane at infinity, and we obtain 



00 



where Resf{uJi) denotes the residuum of f(LU,K) at the pole located at (D^. 

An expression for the propagator in Minkowskian spacetime can be calculated as analytic 
continuation of the euclidean propagator to real times. In this concrete case we insert 

fi-,K) = ^^^l^ (A.26) 

into the frequency summation formula, Eq. (jA.25p . Since r is restrained to the finite interval 
[0,(3], the integration contours can be closed appropriate to Eq. (|A.25p in the complex cj-plane. 
Accordingly, we get 

(2^)4 k^+m^ ^ J {27rr 2E ''^^ '^^ + > 

/j4j, „ifc4T+jk(x-y) C rl^h 1 
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where E = y/kP + rn^ , uj = k4, and nsiPE) = (e'^'^ — 1) ^ denotes the Bose-Einstein distribution. 
We rotate back to real time values by setting r it, — > ik^ and arrive at 

(A.28) 

Inserting 

5(/fc2 _ = ^ {5{k° -E)+ Sik"" + E)) (A.29) 
we arrive at the final result for the propagator of a scalar field in Minkowskian spacetime0 

Ail.(x~y) ^ /(gjAil.(^)e--(^-), (A.32) 
A^L(fc) = k^-L,2+,, +'^^^B{m5{k'~m^). (A.33) 

We have to emphasize that the derivation of the propagator above is merely valid in a free theory. 
In an interacting theory the finite-temperature real-time Feynman rules are more complicated and 
a doubling of degrees of freedom is required to obtain the right perturbative expansion. The 
corresponding propagator is then given in form of a 2 x 2-matrix. 



A. 3 Feynman rules in real-time finite temperature field the- 
ory 

In this subsection we state the Feynman rules for quantum field theory at finite temperature. 

In unitary Coulomb gauge, (j) is diagonal and the pure-gauge ground state field a^^ vanishes 
as mentioned in a previous chapter (unitarity). This physical gauge fixing does not require the 
introduction of Faddeev-Popov ghost fields. Subsequently, the remaining gauge freedom provides 
transversality of the massless gauge mode diSaf^^^ — (Coulomb). A calculation in the real- 
time formalism of thermal field theory is advantageous, since the propagator splits up into a pure 
thermal and a pure vacuum part. In the following we list the Feynman rules: 

The tree-level propagator for a free TLM mode is given as 

P) = ^Sat (Pj. + 2n6ip')ns{P\po\?j - , (A.34) 

where nsix) = ^rzi denotes the Bose-Einstein distribution function and the projector 

P^{p) = <5,,~^. (A.35) 
The tree-level propagator for a free TLH mode is given as 

Dl^,ab{P,P) = -SabD^Ap) (- '-T^ + 2^'5(P' - m^)nBmPo\)) , (A.36) 

' \p^ — + le ) 



^For completeness we mention that in case of a fermionic particle the propagator is given as 

s'iL.(--y) = /(g^41U(P)e--'^-^ (A.30) 

•^mL(p) = ^— -H27rnf,(/3i?)(|* + m)5(p2 (A.31) 

f — ra-'rie 

where np(x) denotes the Fermi-Dirac distribution function np = ^J'^-^ 
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Figure A. 2: The four-boson-vertex and the three-boson-vertex 



with the tensor structure 

bM^(g,,~^-^) . (A.37) 

The three-gauge-boson- vertex and the four-gauge-boson-vertex are the conventional ones 

rjl^, = e{2n)%p + q + k)U,[g^^'^{q-pY+g^P{k-qr+gP^^{p-kr], (A.38) 



[3] ahc 

+facefMe{9^''9''" ~ Q^^g"') + I ade f Me^Q^ " " 9^" 9^] ■ (A.39) 



The four-momenta and cofor indices are defined in figure (|A.2[) . Finally, one has to divide every 
loop diagram by i and the number of its vertices |42j . 

A. 4 One-loop polarization tensor of the massless gauge mode 

In this section we consider the one-loop polarization tensor for the massless gauge mode in SU{2) 
Yang-Mills thermodynamics. This object was first calculated in [IB] and led to some astonishing 
implications [45] . It is evident, that if than elementary particle propagates in a medium, its 
properties are modified caused by interaction with this medium. As we will see, these effects 
dynamically generate an effective mass. In case of SU (2) Yang-Mills thermodynamics this medium 
is given as the nontrivial vacuum, the propagating gauge boson is the TLM-mode 7, interacting 
with the TLH- modes . AU these non-Abelian effects are summarized in the polarization tensor 

The polarization tensor of the TLM-mode is transversal for any value of p^, 
P^n^" = . (A.40) 
Thus the following decomposition in a transversal and a longitudinal part is valid: 

n^"- = G(po,p)pr +^(po,p)pr , (A.41) 

where the projection operators are related via 

P^'' = ^ - g^'' - P^'' , (A.42) 

and P^" defined in Eq. ^EM- 
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Figure A. 3: Diagrams contributing to the one-loop polarization tensor of the massless mode. 



In cuclidcan signature the full propagator of the massless gauge mode results as 



p. 



1 



2^<5(/)ns(/3bo|) 



1 



2 "'/^ '-^y 



(A.43) 



If photon propagation would not be affected by interactions with the virtual massive gauge modes, 
the screening functions F and G would vanish and the propagator (|A.43p transforms into Eq. (|A.34p . 

According to j4j the 00-component of the polarization tensor I{^^{p) in the limit po = 0, p ^ 
corresponds to an electric screening mass (also called Debye mass) rriei- It holds 



F{po = 0, p ^ 0) = -noo(po = 0, p -> 0) 



mil ■ 



(A.44) 



This quantity is diverging in this limit, what implies that static electric fields are completely 
screened by calorons. On the other side it is obvious that in the same limit G{po = 0, p — > 0) 
disappears, so that static magnetic fields are not screened. Furthermore, we realize that on-shell 
noo(p^ ~* 0) remains finite, due to Bose-Einstein suppression, with the consequence that the 
screening function F(ji^ 0) vanishes in that limit. 

Relevant for this thesis is the screening function G{po,p), since G modifies the dispersion law 
for the TLM-mode. Instead of 



(A.45) 



a;2(p) 



p2 + $neG(a;(p),p,A,T) 



7(P) - -— JmG(c^(p),p,A,r) 



(A.46) 
(A.47) 



holds. In this context we can identify the square root of the real part of G as a dynamical mass. 

The polarization tensor 11'^'' was computed in 45J and [46 in on-shell approximation ^ 0. 
We should shortly render the essentials of the calculation as well as the results. 

The one-loop polarization tensor is evaluated as the sum of the two Feynman diagrams depicted 
in Fig. lA.3l An important fact is that for thermal massless modes this on-shell polarization tensor 
is exactly calculated in one-loop order. This is resulting from j471 158] . In this work it was proved 
that the irreducible three-loop contribution to the pressure vanishes identically. If we cut one of 
the massless lines in this three-loop diagram we obtain the irreducible two-loop contribution to the 
on-shell polarization tensor of the massless mode, see Fig. lA.4l Thus the two-loop contribution is 
exactly nil, due to the vanishing three-loop diagram. 
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Figure A. 4: Cutting of one of the massless lines in the three-loop diagram for the pressure yields 
the irreducible two-loop contribution for the on-shell polarization tensor. 



At first we have to emphasize that exclusively the thermal part of the propagator contributes to 
the polarization tensor, each vacuum part is exactly nil due to the constraint on the maximal off- 
shellness of the massive gauge mode, Eq. (j2.96p . in combination with the fact that in the deconfining 
phase e > ^/Sn holds. Furthermore, the calculation shows that the contribution of diagram A is 
identically zero in the on-shell limit by virtue of momentum conservation. 

The tensor structure of the tadpole diagram B evaluates to [H] 



3- 



2TT6{e ~m^)nB{l3\kQ\), 



(A. 



where the four-vertex constraint (|2.97|1 implies severe restrictions on the limits of integration. 
Without constraints on generality, p is assumed to point into the three-direction. Then the diagonal 
components flu and 1122 are equal to G. Momenta and temperature are rescaled to dimensionless 
variables subject to 



T ' 



k 



X = 



2'kT 



Introducing cylindrical coordinates for y, 
the screening function G is calculated as 



(A.49) 



(A.50) 
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(A.51) 



where the integration limits are given as 

Pm = 



7r\2(4e2-l)2 27r4e2-l 
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(A.52) 
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The screening function G was calculated numerically in [45j and is drawn in Fig. IA.5l for various 
temperatures. 

As a consequence of the modified dispersion law (|A.46p photon propagation is affected by 
screening effects arising from the polarization tensor. This effects manifests for temperatures a 
few times the critical temperature. For momenta to the right of the sign change in the screening 
function (associated with the dips in the curves) G is negative and possesses a small modulus. 
Thus the energy of the propagating TLM-mode is reduced in comparison to the free case. This 
effect is known as antiscreening However, for momenta to the left of the sign change, we enter 
the regime of screening, G is positive and reaches sizeable values ~ |pp. The photon acquires a 
screening mass. Thus the photon can travel through the plasma only up to a distance G~^/^. The 
situation is summarized in Fig. lA.5l 

Subsequently, fascinating implications concerning the black body spectrum are arising from the 
modified photon propagation properties. Photon propagation subject to an U{1) gauge symmetry 
leads to the well-known Planck law for the energy density of the black body radiation 



If photon propagation is considered in the framework of SU (2)cmb the energy density modifies to 



where uji and uj2 denote the solutions of the equation cj^ — G{lu) = 0. Photons of frequency 
uJi < Lo < UJ2 possess an imaginary momentum modulus and can not propagate. This regime of 
strong screening causes a gap in the black body spectrum at low temperature |48j , depicted in 



^Nevertheless the influence of the tree-level massive modes V is very weak due to the strong restrictions onto 
the maximal ofT-shellness. 




(A.53) 




(A.54) 



Fig.Ea 
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Figure A. 5: log^Q |^| as a function of |p| for A = 1.12 Ac (black), A = 2 Ac (dark grey), A = 3Ac 
(grey), A = 4Ac (light grey), A = 20 Ac (very light grey). This result is obtained by appealing to 
the approximation — 0. The full calculation shows similar results for finite |p|. However, there 
we have lim|p|_,Q |^| > in contrast to the here-indicated result. The dashed curve is a plot 
of the function /(|p|) = 21og]^Q |p|. Here A = 13.87^ = Photons are strongly screened at 

|p|-values for which logj^Q |^| > /(|p|), that is, to the left of the dashed line. The dips correspond 
to the zeros of G. 
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Appendix B 



Technical details concerning 
Chapter 3 



B.l Wick decomposition of the correlator 

This section contains the Wick decomposition of {&oo{x)Qoo{y)) into a sum of products two-point 
functions of the elementary gauge fields Afj_{x). The evaluation is not a mathematical challenge, 
but nevertheless a cumbersome task. We emphasize that our focus is on the connected correlations, 

local diagrams are neglected. 



(eoo(a;)eoo(2/)) = 

= {{-Fox{x)Fo \x) + ifF^,{x)F^\x)){^For{y)Fo ^{y) + 3-fF,r{y)F^%y))) 
= {Fox{x)Fo\x)For{y)Fo ^y)) - if{Fox{x)Fo \x)F,r{y)F-^ {y))) 

^^{F,,{x)F^\x))For{y)Fo^iy)) + ^^{F,x{x)F^Hx))F,r{y)F^^y))) 
= {{d,oA\x)){d,oAx{x)){dyoA^y)){dyoAr{y))) - {{d,oA\x)){d,oAx{x)){dyoA^{y)){dy.Ao{y))) 
-{{d,oA\x)){d,oAx{x))idy^Ao{y))idyoAr{y))) + {{d,oA\x)){d,oAx{x))^^^^ 
-{{d,oA\x))id,.Ao{x))idyoA-{y))idyoA,{y))) + {{d,oA\x)){d,.A„ix)^^^^ 
+ {{d,oA\x))id,.Ao{x)){dy^Ao{y)){dyoAr{y)))-{{d,oA^{x))id,.Ao{x)^^^^ 
^ ( (a,, Ao (x) ) (9,0 A;, (a;) ) (a,o (y ) ) (a,o A. (y ) ) ) + ( (9,, Ao (.t) ) (9,0 A;, (a;) ) (a,o A- (y ) ) (a,. Ao (y ) ) ) 
+ ( (a,, (.t) ) (9,0 A;, (x) ) (a,, (y ) ) (9,0 A. (y ) )) - ( (a,, (.t) ) (9,0 A;, (x) ) (a,, (y ) ) (a,. Ao (y ) )) 
+ ((a,,Ao(x))(9,xAo(x))(a,oA-(y))(a,oA,(y)))-((9„Ao(x))(9,xAo(x))(a,oA-(y))(a,.Ao(y))) 



X [( (9,0 A\x) ) (9,0 (x) ) (9,„ (y ) ) (9,,. (y ) ) ) - ( (9,0 (x) ) (9,0 A;, (x) ) (9,„ A- (y ) ) (9,. A„ (y ) )) 
- ( (9,0 (x) ) (9,0 A;, (x) ) (9,, A'^ (y ) ) (9,. A. (y ) )) + ( (9,o (x) ) (9,o A;, (x) ) (9,, A^y)) (9,. A,{y))) 
-((9,oA\x))(9,xAo(x))(9,,A-(y))(9j,.A,(y))) + ((9,oA\x))(9,.Ao(x))(9,<,A-(y))(9,.A<,(y))) 
+((9,oA\x))(9,.Ao(a;))(9,,A-(y))(9,.A.(y)))-((9,oA^(x))(9,.Ao(x))(a^,A'^(y))(9,.A<,(y))) 

-((9,,Ao(x))(9,oA;,(x))(9,„A,(y))(9,.A,(y))) + ((9„Ao(x))(9,oA;,(x))(9,„A,(y))(9,.A,(y))) 
+ ( (9„ (x) ) (9,0 (x) ) (9,^ A, (y ) ) (9,. A, (y ) )) - ( (9„ Ao (x) ) (9,0 (x) ) (9,^ A, (y ) ) (9,. A, (y ) )) 
+ ( (9,, Ao (x) ) (9,. Ao (x) ) (9,„ A. (y ) ) (9,. A. (y ) )) - ( (9„ Ao (x) ) (9,. Ao (x) ) (9,„ A. (y ) ) (9,. A. (y ) )) 
-((9„Ao(x))(9,.Ao(x))(9,,A.(y))(9,.A,(y))) + ((9„Ao(x))(9,.Ao(x))(9,,A,(y))(9^.A<,(y)))] 
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( 



-.900 

4 



x[{{d,^A\x)){d,.Ax{x)){dyoA^y)){dyoAriy)))-{{d,^A\x))id,.Ax{^^^^ 
-{{d,^A\x)){d,.Ax{x))idy^Ao{y))idyoAr{y))) + {{d,^A\x)){d,.Ax{^^^^^ 
-{{d,^A\x)){d,.A^{x)){dyoA^{y)){dyoAr{y))) + {{d,„A\x)){d,.A^{x^^^^ 

M{d.A\md.>^A,{x)){dy^Ao{y)){dyoAAy)))-{{d,^A\x)){d,.A^{x))^^^^ 
-{{d,,A^{x)){d.,.A^{x)){dyoA^{y)){dyoAr{y))) + {{d..,A^i^^^^^ 
Hid..A-{x)m^A,{x))idy^Ao{y))idyoAriy)))-{(d,,A-{x))id,.^ 
M{d.,A-{x)){d,.A,{x)){dyoA-{y)){dyoA,{y)))-{{d,,A'^{x)){d,.A^^^^^^^^ 
-{{d,,A'^{x))id,.A,{x))idy^Ao{y)){dyoAr{y))) + {{d,,A'^{x)){d,^^^ 



x[{{d,^A\x)){d,.A^{x))idy^A^y)){dy.AAy))) ~ {id..A\x)){d^ 
-{id,^A\x))id,.A^ix)){dy^A'^iy)){dy.Ariy))) + {id,^A\x)^^^^ 
-{id,^A^ix)){d,.A,{xmdy^A^y))idy.Ar{y))) 

+m^A\x)){d,.A,{x))idy^A^{y))idy,Ar{y)))-{{d,^A\x))id,^^^^ 
-{{d,,A'^{x)){d,.A^ix)){dy,A^y)){dy.Ar{y))) + {id,,A^{x)){d,^ 

Mid..A-ix)){d,^A,{x)){dy^A''{y)){dy.AAy)))-{{d.,A-{x)){d,.A,{^^^^^ 

M{d.,A'^{x)m,.A^{x)){dy^A^{y)){dy.AAy)))-{{d.,A-{x)){d,.A,^^^^^ 

-{{d,,A-ix)){d,.A,{xmdy^A'^iymdy.AAy))^ 



{id,oA\x)){d,oAx{x)MdyoA^y))idyoA,{y)))+2{id,oA\x)){dyoA^y)^^^^^ 

-{id,oA\x))id,oAxix))){idyoA^y))idyMy)))-M{d.oA\x)){dyoA^y))){^^^^^ 

-{id^oA\x)){d,oAx{x))){idy^Aoiy))idyoAr{y))) 

Mid.oA^{x)){d,oA^{x))){{dy^Ao{y)){dy.Ao{y)))+2{{d,oA^{x)){dy^A 
-{{d,oA\x)){d,.Ao{x))){{dyoA^y)){dyoAr{y)))-2{{d,oA^{x)){dyoA^y)^^ 
+ {id,oA\x)){d,.Aoix))){idyoA^y))idy.Aoiy)))+2{id,oA^ix))(dyoAr(^^^^ 

M{d.oA\x)){d,.Aoix))}{idy^Ao{y)){dyoAr{y)))+2{{d,oA\x)){dy^Ao{ym 
-{{d,oA\x)){d,.Ao{x)Mdy^Ao{y)){dy.Ao{y)))-2{{d,oA\x)){dy^Ao{^^^^^ 
-{{d,,Ao{x)){d,oAx{x))){{dyoA^y)){dyoAr{y))-2{{d,,Ao{x)){dyoA^y))){{^^^ 

+ {{d^Mx)){d,oA^{x))){{dyaA^{y))idy.Ao{y))+2{{d,Mx)){dy^Ao{^^^^ 

+ ( (a,, Ao (x) ) (9,0 A A (a;) )) ( (a^^ Ao (y ) ) (9,0 A, (?y)) + 2 ( (a,, Ao (x) ) (9,0 (y ) ) ) ( (9,0 A A ) (a,. Ao (?y) )) 
-((9,,Ao(x))(9,oAA(x)))((9,,Ao(y))(9,.Ao(?y))-2((9,,Ao(x))(9,,Ao(y)))((9,oAA(x))(9,.Ao(?y))) 
+ ((9,,Ao(x))(9,xAo(x)))((9^oA^(y))(9,oA,(y))+2((9,,Ao(x))(9,oA^(y)))((9,x^o(a;))(9yoA,(y))) 

-((9,,^o(a;))(9,.Ao(a;)))((9,o^-(y))(9^.^o(2/)))-4((9,,^o(a;))(9,oA-(j/)))((9,.Ao(x))(9^.^o(y))) 

-{{d,,Ao{x)){d,^Ao{x))){{dy^Ao{y)){dyoAriy))) 

+((9,,^o(x))(9,.Ao(a;)))((9,,Ao(y))(9^.^o(y))) + 2((9,,Ao(x))(9^,^o(y)))((a,A^o(a;))(9,.Ao(y))) 




(B.l) 
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X [{{d,.A\x)){d,oA^[x))){{dy^A^{y))[dy.A,{y)))+2{{d,oA\xm^^^ 
-{{d^oA\x)){d^oA^{x))){(dy^A^{y)){dy.A,{y)))-^{d^oA\x)){dy^A% 
- ((9,0 (x) ) (9,0 A^ (x))) {{dy^ A'^ (y) ) {dy. A, (y))) 
Mid,oA\x)){d,oA^{x))){{dy^A'^{y))idy.A,{y)))+2{{d,oA\ 
-{{d,oA\x)){d,.Ao{x))){{dy^A^y)){dy.Ar{y)))-2{{d,oA\x)){dy^A-{ym^^^ 
+((9,oA^(a;))(9,.^o(a;)))((9,,^-(y))(9^.^.(2/))) + 2((9,o^^(a;))(9,,A'^(y)))((9,.Ao(a;))(9,.^,(y))) 
+ ((9,oA^(x))(9,.Ao(a;)))((a,,A'^(y))(a,.^,(y)))+2((a,oA^(x))(a,„A-(?y)))((a,^Ao(x))(9,.A^(y))) 
-{{d^,,,oA\x)){d,.Aoix))){{dy^A'^{y)){dy.A^{y)))-2{{d,oA^^^^^^^ 

-((a,,Ao(x))(a,oA,(x)))((a,„A-(y))(a,.A.(j/)))-2((a,,^o(x))(a,^A-(y)))((a,oAA(.x))(9,.A,(2/))) 
+((a,,Ao(x))(a,oA,(x)))((a,,A-(y))(a,.A<,(y))) + 2((a,,Ao(x))(a,,A'^(y)))((a,oA,(x))(a,.A,(y)^ 

+((9,,Ao(a;))(9,oAA(x)))((a,,^-(j/))(a,.A,(y)))+2((5,,Ao(x))(9,<.A-(2/)))((a,oA;,(a;))(a,.^,(^ 
-{id,,Aoix)){d,oAx{x))){idy^A^y)){dy.A^{y)))-2{{d,,Ao{x)){dy^A^^ 

+((a,,Ao(x))(9,.Ao(x)))((a,„A-(y))(9,.A,(y)))+2((a,,Ao(x))(a,.A-(j/)))((a,.Ao(x))(5,.A.(t/))) 
^((a,,Ao(x))(a,.Ao(x)))((a,„yl.(j/))(9,.A,(y)))-4((5,,Ao(x))(a,^^-(y)))((a,.Ao(x))(a,.A,(2/))) 
^ ((a,, Ao (x) ) (a,. (a:))) ( {dy^ A- {y)){dy. A^ (?;))) 

+ ((a,,Ao(x))(a,^Ao(a;)))((a,,A'^(y))(a,.A.(y)))+2((a,,Ao(a;))(a,^A,(y)))((a,xAo(a^))(5,.A,(y)))] 



X [((5,„^^(x))(9,.AA(x)))((a,o^-(2;))(a,oA,(y)))+2((a,„A^(x))(a,oA-(y))>((a,.AA(x))(a,oA,(y))) 

-((a,,A^(x))(a,»A;,(x)))((a,oA-(y))(a,.Ao(y)))-4((a,.A^(x))(a,,Ao(y)))((a,.A;,(x))(a,oA,(y))) 

-{{d,^A\x)){d,.A;,{x))){{dy^Ao{y)){dyoAriy))) 

+((a,„A\a.))(a..^A(a;)))((a,^^o(y))(a,.^o(2/)))+2((a,,A^(a;))(a,^Ao(2/)))((a,.A,(a.))(9,.^o(y))) 
-((a,„A\a;))(a,.A„(a;)))((a,oA-(2/))(a,o^,(y)))-2((a,^A\a;))(a^o^-(2/)))((a,.A«(a;))(aj,o^,(y))) 

+((a,,A^(a;))(9,.A«(a;)))((a,oA-(y))(a,.Ao(y)))+2((a,„A\a;))(9^,Ao(y)))((a,.A„(x))(a,oA,(y))) 

+((a,,A^(x))(a,.A«(x)))((a,,Ao(j/))(a,oA.(y))) + 2((a,.A^(x))(a,oA-(y)))((a,.A«(x))(a,.Ao(y))) 
-((a,„A^(x))(a,.A„(x)))((a,,Ao(y))(a,.Ao(y)))-2((a,,A^(x))(a,,Ao(?y)))((a,.^«(.T))(a,.Ao(y))) 
-((a,,A«(x))(a,.AA(x)))((a,oA-(y))(a,oA.(y)))-2((a,,A«(x))(a,oA-(y)))((a,.^,(x))(a,oA.(y))) 
+((a,,^«(x))(a,.A;,(x)))((a,oA-(y))(a,.Ao(y))) + 2((a,,A«(x))(a,,Ao(y)))((a,oA;,(x))(a,oA,(y))) 
+((a.,^«(x))(a..^A(a;)))((a,^^o(y))(9,oA.(2/))) + 2((a,,^«(x))(a,a^-(y)))((a.o^A(a;))(a,.Ao(2/))) 
-((a,,^'^(x))(9..^A(a;)))((a,^Ao(y))(9,.^o(2/)))-2((a,,A«(x))(a,^Ao(2/)))((a,oA,(x))(a,.Ao(y))) 
+((a,,A«(x))(a,.A.(x)))((a,oA-(y))(a,oA.(y))) + 2((a,,A''(x))(a,oA-(2;)))((a,xA«(x))(a,oA.(y))) 
-((a,,A«(x))(a,.A.(x)))((a,oA-(?y))(a,.^o(?y)))-4((a,,^'=(x))(a^o^-(t/)))((a,.^„(x))(a,.Ao(y))) 
- ((a,, (x) ) (a,. A« (x))) ((a,, Ao (y) ) (a,o a, (y) )) 

+((a,,A'=(x))(a,.A.(x)))((a,,Ao(y))(a,.Ao(y))) + 2((a,,A«(x))(a,,Ao(y)))((a,.A,(x))(a^.Ao(y)))] 



X [((a,^^^(x))(a,.^A(x)))((a,„A-(y))(a,.A,(y))) + 2((a,,A^(x))(a,^A-(y)))((a,.AA(x))(a,.A.(j/))) 
-((a,„A^(x))(a,.AA(x)))((a,„A-(y))(a,.A.(2/)))-4((a,„^^(x))(a,^A-(y)))((a,.^A(a;))(a,.^,(y))) 

-{{d,^A\xym.A^{x))){{dy^A^{y)){dy.Ar{y))) 
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local diagrams real correlations 



(B.2) 

We are interested in the connected two-point correlation function, thus we neglect the local 
diagrams in Eq. (|B.2[) . Inserting the propagator for the massless gauge mode in Coulomb gauge 
and summarizing the terms of the same structure leads us to Eq. p.Sp : 

{&oo{x)eoo{y)) - 2{d,aAAx)dyoAAy)){dxoAAx)dyoAAy)) 

-goo (9,0 A^ (x)dy. A^ {y)){d,o A^ {x)dy^ A^ {y)) 
+500 {dxo A^ {x)dy, A^ iy)){dxO Ax {x)dy. A, (y)) 
-goo {dx^ A^ {x)dyO A^ (y)) {d^. Ax {x)dya (y)) 
+500 {dx^ A^ {x)dyO A^ {y)) (5,a A^ {x)dy. AAy)) 

+ ^ (9,^ A^ {x)dy, A^ (y)) {d.,. Ax {x)dy. A^ (y)) 
- ^ (dx.A^ {x)dy^ AAy)) {dx"Ax {x)dy. A, (y)) 
-^{dx^ A^ {x)dy^ A^ (y)) (9,. A^ {x)dy. A, (y)) 

+ ^ {dx^A^ {x)dy^ AAy)) (5,A A, {x)dy. AAy)) 

+2{dxoAAx)dyoAAy)){dxoAoix)dyoAo{y)) 
-2{d,oAAx)dy.AAy)) (5,0 Ao {x)dy^Ao{y)) 
+2(5,. A° ix)dy. A° (y )) (9,^ Ao {x)dy^ Ao (y)) 
-2goo{dx^>AAx)dyoAAy)){dxoAo{x)dyoAo{y)) 
+4500 (9,0 AO {x)dy. AO (y)) (9,0 Ao {x)dy^ Ao{y)) 
-2goo{dx^AAx)dy.AAy)){dx^Ao{x)dy^Ao{y)) 

+ ^ (9,0 AO (x)9,o AO (y )) (9,0 Ao (x) 9,0 Ao (y )) 

-5oo (9.0 AO {x)dy. AO (y)) (9,o Aq {x)dy^ Ao (y )) 

+ f^(9,.AO(x)9,.AO(y))(a,^Ao(x)9,„Ao(y)). 

(B.3) 

Subsequently, we take the partial derivatives out of the expectation value and insert the propa- 
gator for the massless mode, Eq. (IA.34p . 
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B.2 Calculation of the thermal contribution in Coulomb 



gauge 

B.2.1 The entire thermal part including oscillatory terms 

The evaluation of the integrals over the zero-components of the momenta eliminates the Dirac 
^-distribution and yields additional oscillatory terms. Since we are focussed on a time-averaged 
expression for the energy transport, we are setting — if for the subsequent calculation. At this 
point we state the time-dependent two-point correlator after performing the integration over po 
and k[j for reasons of completeness. Starting with Eq. I|3.6p we insert the decomposition Eq. (|3.12p 
and perform the integration by employing the convolution property Eq. (|3.1ip of the (5-distribution. 



(eoo(a^)eoo(tj))*'' 



■{5{po - IpI) + Sipo + \p\)mko - |k|) + Sih + |k|)) 



iS{po - IpI) + -5^0 + |p|))(<5(fco - |k|) + S{ko + |k|)) 
.00 / dpo J dko J / ( (^) ' - l) pofcopk . nM) n^kol) • e--(^-)e-(^-) 



■Wpo - IpI) + -5(^0 + |p|))('5(fco - |k|) + 5(fco + |k|)) 



iSipo ~ IpI) + s{po + |p|))(<5(fco - |k|) + s{ko + |k|)) 

+ .00 / dpo J dko J J ( (^) ' - l) .ofcopk . nM) n,iP\ko\) • e--(^-)e^-(^-) 

■{Sipo ~ IpI) + s{po + |p|))(<5(fco - |k|) + s{ko + |k|)) 

+ f /..o / dk. I ^ / (l + (^) ipkr • nmP.\) n.imol) . e-^^(^-)e-(^-) 

■{S{po - |p|) + S{po + \p\)mko - |k|) + S{ko + |k|)) 
- .0^0 / dpo J dko J / ( (^) ' - l) Pkp. . uM) nsiPlkol) . e-(^-)e--(^-) 

■{Sipo - IpI) + 6{po + |p|))(<5(fco - |k|) + S{ko + |k|)) 

. f >->^o/(^/(^ ((^)^l)^^k^...^ 
•(<5(po - IpI) + -^(Po + |p|))('^(A;o - |k|) + Sik, + |k|)) 

(B.4) 
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(27r)3 J (27r)3 \ VlPl|k 

•0 

d^p f d^k I / pk ^ ^ 



j |p||k| • nB(/3|p|) nsim) ■ e^''^""''^ e^'^^''"''^ 
(cos((|p| + |k|)(x - v)^) + cos((|p| - |k|)(x - y)o)) 
p||k|.nB(/3|p|)ni3(/3|k|).e'P(^"y) e^'^'^-y) 



(27r)3 7 (27r)3 \ VlPl|k| 

•(cos((|p| + |k|)(x - v)n) + cos((|p| - |k|)(a; - y)o)) 
d^p /■ d^k I / pk \ / pk ^ ^ 



|p||k| • nB(/3|p|) nB(/3|k|) • e^P'^'-y) e^'^^''-^) 



(2^)37 (27r)3 \^VlPl|k|; VlPl|k| 

•(cos((|p| - \k\){x - y)o) - cos((|p| + \k\)ix - y)o)) 
d^p f d^k I / pk \^ f pk 



\p\\k\ ■ nB{/3\p\)nB{/3\k\) ' e'^^""'') e^'^^''-^) 



(2^)3 y (27r)3 \^VlPl|k|; VlPl|k| 

•(cos((|p| + \k\){x - y)o) ~ cos((|p| - |k|)(a; - y),)) 
d^p f d^k ( / pk \^ / pk 



47 (2.)3 J (2.)3 - viPiiki; ' viPiiki; ; ipiw--(/^ipi)-(/'W)-e-^^-''^e-(-) 

•(cos((|p| + \k\){x - y)o) + cos((|p| - |k|)(a; - y)o)) 

•(cos((|p| - \k\)ix - y)o) - cos((|p| + \k\)ix - y)o)) 

•(cos((|p| + |k|)(x - v)^) + cos((|p| - |k|)(a; - j/)o)) 

•(cos((|p| + \k\){x - y)o) - cos((|p| - |k|)(a; - y)o)) 
+ (^)') IPl|k|-n.(^|p|)n.(/^|k|).e^^(-^)e-(-^) 

•(cos((|p| + \k\){x ~ y)o) + cos((|p| - |k|)(.T - y)o)) . 

(B.5) 

Restricting to equal-time and summarizing the terms of identical structure leads us to Eq. p.l3p 
(eoo(x)eoo(y))"' = (y'^|p|nB(/3|p| 

+ (^)^IPlW-B(^|p|W(/3|k|)e-e-, 



2 

I e'P^ ' 

2 



(B.6) 



B.2.2 Integration over azimuthal and polar angle 



In this subsection concerning the thermal contribution we proceed with the detailed calculation 
of the integrals over the azimuthal angle (j) and the polar angle 6. Introducing three-dimensional 
spherical coordinates for the momenta p and k we formulate the spatial scalar product pk in terms 
of the azimuthal and polar angles. We evaluate the integration over the azimuthal angles cjjp and 
(j}k- As a consequence the integrals factorize in the remaining variables radial momentum-modulus 
and polar angle. We obtain 
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(eoo(x)eoo(y))"' 



oo +1 \ ^ 



\0 -1 

C30 +1 OO +1 



+2Tr- Jd\p\ JdcosepJd\k\ JdcosOk J^^^_ 1^^^ ^ ei|p||z| cos^^ gi|k||z|cose. 



-1 

2tt 



(elPl - 1) (elk| - 1) 
d(j)p (sin 9p sin cos (pp + cos cos 9k)^ 





oo +1 



47r2( ld\p\ Idcose 'P'' .'IPlNcosg 



(27r)6/38 I y y (elPl - 1) 

\ \0 -1 

y^lpl cos 0sin2e^^^L_e»lpl|S| cose 



vO -1 

oo +1 



W I y^lpl ydcos^cos2 0^^^^e'IPll^l=°^^ 

-1 



(B.7) 



Subsequently, we expand the exponential into spherical harmonics and express the powers of 
COS0 in terms of the first Legendre polynomials Pi(cos0). The Legendre polynomials up to the 
fourth order are given as: 

Po(cos6i) = 1, 
Pi (cos 6) = cos , 

P2(cos6l) = ^(3cos^6l- 1), 

Pg (cos 6») = i (5 cos^ 61 - 3 cos 6*) , 

P4(cos(9) = ^(SScos^fii-SOcos^ei-S). 

(B.8) 

Employing their orthogonality relation 
+1 

jdxPm{x)Pn{x) = l^^^^^mn ■ (B.9) 
-1 

leads us to 
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(27 



oo +1 



(eoo(x)eoo(y))' 

(oo +1 ^ 

j d|p| y" dcos^-|^^^'(2/ + l)J^(|z||p|)P^(cos0)Po(cos( 
-1 

? I ~ |3 °° 

d|p| / dcos^^^L^^z'(2; + l)ji(|z||p|)P/(cos0)Po(cos0) 

-I Jd\p\ J dcos0^^^Y.t'i2l + lln{\i\\p\)Pi{cose)P2{cose) 
-1 '=0 y 

(oo +1 ^ 

i y d|p| J dcos0^^^^z'(2; + l)jK|z||p|)P/(cos0)Po(cos( 
-1 

V I ~ |3 

d|p| / (icos0^^^^z'(2; + l)jz(|z||p|)P/(cos0)Po(cos( 

^|Pld^^°(|Pll^l) + ^/ '^|Pl^[^J2(|p||i|) 



27r2. 



• I o / d|p|^f^Jo(|p|N)-- / d|p|^^j2(|p||z|) 



IPP 



3iPi _ r 



(B.IO) 



Summarizing integrals of the same structure let us finally arrive at Eq. p.2ip 

2 



(900 (x) 000 (y))"' = 



1 / 647r^ 



(27r)6/38 \ 3 
32^2 



rf|p|^|r^j2(|p||i|) 



(B.ll) 
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B.2.3 Integration over momentum 

The remaining integrals involve spherical Bessel functions of the first kind. These function are 
related to the ordinary first kind Bessel functions of fractional order via 



The (spherical) Bessel functions satisfy the following recurrence formula: 
The spherical Bessel functions for = 1 , . . . , 4 are given as 



3q{z) = 


sm z 

z 


Jiiz) = 


sin z 




(-- 

\z3 


hiz) = 


(S- 


jiiz) = 


/105 



105 45 1 



15 1 



cos z , 
/105 10 



COS z . 



In case of the thermal contribution we employ the two integral formulas [44] 



dx 



(-1)^ 



a2m 







|coth(7r6)-l 



„2m+l 



dx- 



cos{bx) 



Q2m+l 







where 6 > 0, m > and ni G Z, and obtain 



|coth(7r6)-l 



(B.12) 



(B.13) 



(B.14) 

(B.15) 
(B.16) 



elpl - 1 



^/^IP 



1 

jii4 



^, |p psin(|p||i|) 
elPl - 1 

TT'^ coth(7r|z|)cosech^(7r|z|) 



(B.17) 



''IPI elPl - 1 



^1 |Sinam_ 1 
elpl-l z 



.,|Ppsin(|p||z| 



elpl - 1 











d|p| 



IpI cos(|p||z|) 



cIpI - 1 

-8 + 7r|z|(3coth(7r|z|) + 7r|z|(3 + 27r|z| coth(7r|z|)) cosech^(7r|z|)) 







2|z 
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(B.18) 
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B.3 Calculation of the vacuum contribution in Coulomb 



gauge 

B.3.1 Integration over azimuthal and second polar angle 

As similar procedure is applied to the case of the vacuum propagator. Again Eq. (|3.5p provides a 
benchmark for the calculation. In contrast to the thermal contribution we employ a Wick-rotation 
to euclidean signature (po7^o ^ —ipo, ~iko, x^,y'^ —ix'^,—iy^ and g^^i, —Sf^u). Thus we 
are able to transform the variables to 4D spherical coordinates and choose to point into the 
3-direction. In detail we have 



(27r)4 J (27r)4 ''^a 1,2 



(27r)4 J (27r)4 p2 ^2 



(2^)4 J (2^)4^ 

^™ y (2^)4 y (2^)4^ ^^'^ ^2^2^ e 
5|, /• d^p f d^p pXTp^^ P^P" kxkr ^^^pc^-ikC 



2 J (27r)4 y (27r)4 "° p2 



(27r)4 7 (2^)4^^"-" ' ^ p2 k2 

M^P . 2,2 o , ,r , ,A ,.,e-^PCe-^'=« 



i2 /■ ^4„ /• ^ p-ipC p-ifcC 



rf-^p f d^p ^ ^ pop„ kokr _,p(; 



k2 



2 7 (27r)4 j {2nf p^ k^ 

"° J (27r)4 y (2^)4 p2 ^ 



2 y (27r)4 y (27r)4 p2 ^2 

(2^ J -(^(P^^'^O-PokoPrkr+PxkxPrkr)-^-^ 

/d^p f d'^p e*P'» e'*"» 

(2^ y (2^'-^°^° ~ '^PokoPrkr + PxkxPrkr) 

+ 1"^ y (2^ y (^(?'ofco-2pofcoPrfcr+PAfcAPrfc.)-^-j^ 



(B.19) 
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2U Wr^'^j w'^' ^ 



+6 



d'^p f d'^k f pk 



(27r)4 y (27r)4 V|p|!k| 
1 f d^p f d'^k f pk 



Pofco|p||k| 



p2 ^2 



2 y (27r)4 y (27r)4 V|p||k| 



+2 



+2 



d^p /■ (i'^fc / pk 



(27r)4y (27r)4 VIp||k| 
(i'^p f d'^k f pk 



Pofco|p||k 



2 J (27r)4 y (27r)4 V|p||k| 



(B.20) 

Introduction of dimensionless variables via rescaling of the momenta and transformation to hyper- 
spherical coordinates subject to the convention 



p= \p\ 



/ cos ip \ 
sin sin cos (p 
sin tp sin sin (f) 
\ sin cos 9 J 



with 



< ■0,6' < TT 

< < 27r 



(B.21) 



|p| = Vaf+xf+af+ar 



yields 

(eoo(x)eoo(y))' 



^^^^1 2 -4^' ( /d|p|dVpd^^p|p|'sin2V;pC0s2^pSin(^pe'l^llfl^'"^'"=°^^'' 



+ 



^ • 277 y <i|p|dV'pC?0p y d|A:|c?V'feC^^fc sin^ Vp cos^ Vp sin 9p\kf sin^ Vfc cos^ tpk sin 



^ilpllfl sin \jjp cos6pgi|fe||C| sin-i/iii cosS^ 
f-27r 



■ / (sin^ 9p sin^ ^fe cos^ ?!'p + 2 sin 6p sin ^fe cos 9p cos + cos^ 9p cos^ ^fc) 

+6 • 27r / d\p\dtppd9p / dlfcldV/tc^^fclpI"^ sin^ ^/jj, cos'i/'p sin0j,|fc|^ sin"^ f/j/c cos^/ifc sin^^ 



. g»|p|KI sini/'p cosepgi|fc||C| sini/n, cosS/i 

■ / d(jf)p (sin 6p sin ^fe cos (^p + cos 6p cos ^fc ) 
Jo 



+ — • 27r y d\p\d-)ppd9p J d\k\d-)pkd9k\p\^ sin^ -tppSinOplkl''' sin'' sin 

. g*|p|ICI Bin t/jp cos fpgilfcllCI sini/ifc cosSfc 
/.27r 

• / d(^p (sin^ 9p sin^ cos^ ^iip + 2 sin 6p sin cos 9p cos 0^ + cos^ 9p cos^ ^fe) 
+2 • 47r2 (i|p|#pci6'p|p|3cos2 V'pSin6lpe^lPll^"l^'"'^''™^^''^ 

+2 • 27r J d\p\dippd9p j d|fc|rfV'ferf^fe|p|^sin^pCOsV'pSin^p|fc|^sin^fe cos^'fe sin^fc 



63 



. p^IpIKI sin V'p cosepgi|fc||C| sinj/'ib cos^ib 

dcpp (sin Op sin 9k cos 0), + cos 6p cos 9k ) 



• 27r d\p\dtppd9p J d\^d^pkd9k\'pf sin'^ 'ippSm9p\k\^ sin^ -ipk sinOk 



^i\p\\(\ sintpp cos6p^i\k\\C\ sinj/i/t cosSfc 



/ d(j)p (sin^ 0p sin^ 0^ cos^ 0p + 2 sin sin 9k cos 0p cos 0^ + cos^ 9p cos^ 
Jo 



(B.22) 



Again the expression above factorizes in two identical contributions for each momenta. Thus we 
are able to write 



(eoo(x)eoo(y))^^" = 



1 9 



CO TT 77 



(27r)8/38 2 



-•47r^ d\p\ d-tp d9\p\'^sm-^ijjcos^ijjsm9e 



aMPWCl sin ip cos 6 



^0 



OO TT 7r 



+ 



i-47r2 ljd\p\ Jdy: ^ri^ [pl^sin* Vsin^^e'IPll^>'"^^°'' 



OO 7r 7r 

Jd\p\ jd^ Jc 



+ - ■ 



d9 \p\^ sin^ ip cos^ ip sin sin"* 9 e 



3 flf,i|p||CI sin'0cos6 



VO 

OO TT 7r 



+ 



^•47r2 jyd|p| Jdi; yd6'|p|^sin2Vcos2Vsin6icos2 6ie'l^*llfl^'"*'=°^^ 



VO 

OO 7r 7r 



+6-47r2 d\p\ Idip /d6l|p|3sin3VcosVcos6ie'lpll^'l"'"*™"^ 



VO 

OO 7r 7r 



+ --27r2 I /# /d6'|^|3sin4V'sin3 6'e*IPll«l''''^'^'=°''^ 



VO 

OO 7r TT 



+ 



^•47r2 jyd|p| y# Jd9\pf sm^TPsmOcos'^ 9 e'\p\\^\^'''''"'°^^ 



KO 

OO TT 7r 



+2 



47r2 jyd|p| y# yd6i|p|3cos2Vsin6le'lpll^'l"'"''"^°^^ 

\0 



OO 7r 7r 



+2-47r /(iT/; /c?6' 1^1*^ sin cos sin 6* cos 6*6 



ICI sintp cos t 



VO 

OO TT TT 



+ /'^IpI /'^V^ /d^lpl^sin^Vsin^^e'l^ll^'l^'"^'^"^^ 



vO 

OO TT TT 



+ ^•277^ ly^l y# j d9\p\^ sir? ipsinOcos^ 9 



i\p\\C,\ sin^/jcos^ 



VO 



(B.23) 
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In order to proceed with the 6'-integration we expand the plane wave into spherical harmonics. 
We replace sin^ ^ by (1 — cos^ 6) and — d^sin^ by dcos^. The polynomials in cos^ are written in 
terms of Legendre polynomials. 
We obtain 

/ / n +1 ^ 

IStt^ ijd\p\ jdi)\p\^sm^ipcos^ip j d cos 9 Pq {cos 0) ^ i\2l + 1) ji{\p\\C\ sin 7p) Pi {cos 6) 
y Vo -1 '=0 

OO TT +1 ^ 

+2Tr'^^jd\p\ jdtp\p\^ sin^ y"dcos6'^z'(2; + l)j,(|p||C| sinV')P;(cos6') 



1 

(27r)8/J8 



VO 



(OO 7r +1 

Jd\p\ Jdiplpf sin'^ tpcos^ tp jdcos0 
-1 



^Po(cose) - |P2(cos^) ) ^J{21 + l)ji{\p\\C\ sin V)P/(cos^) 



3 ' '3 

OC' TT +1 



+187r2 ijd\p\ Jdiplpfsin'^ipcos'^tp jdcosO 



\0 



1 2 \ °° 

-Pq{cos6) + Mcose) ^i'(2Z + l)ji{\p\\C\ sin V)P;(cos^) 



3 ^ ' 3 

OO TT +1 

+2471^ I fdlPl /(iVlpl^sin^V'COsV' /"dcos6'Pi(cos^) ^i'(2^ + l)ji(|p||C| sin7/;)P;(cos^) 



VO -1 

(OO 7r +1 

Jd\p\ Jd^\p\^sm'^ip jdcosO 
-1 



;=o 



^Po(cose) - |P2(cos0) ) ^i'(2Z + l)ii(|p||C|sinV)P/(cos^) 



3 ^ '3 

OO TT +1 



+27r^ (/^l^l j di)\p\^ sm^ i) jd cose 

-Po(cos0) + -P2(cos0)j YJi'^^ + IbKI^IICI sinV)Pi(cos( 

OO +1 OO ^ 

+87rMyd|^| jdip\p\^ cos^ip jdcose^i\2l+l)ji{\p\\l\sm'il))Pi{coae) 
\o -1 '=0 J 

+877^ ( /^^l^l Jd^\pf simp cos 2p J d cos 6 Pi {cos 6) {21 + 1) ji{\p\\C\ sin ■tp) Pi {cos 9) 



VO 



VO 
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OO TT +1 



+977^ ijd\p\ jdil)\p\^shi^ip jdcosO 

KO -1 

^Po(cos0) - ^P2(cos0)^ f2i\2l + l)ji{\p\\C\sm,p)Pi{cos0)^^ 
+1 

+187r2 ijd\p\ jd^plpfsin^^ jdcosO 

^Po(cos^) + ^P2(cos^)) Y.i\2l + l)ji{\p\\C\ simlj)Pi{cose) 



^0 



;=o 



(27r)8/38 



OO TT 



OO TT 

727r^ I y"d|p| jdtp\pfsm^tpcos'^tpjo{\p\\C\smi;) 

^0 

+877^1 yd|p| ydVIPI'sin^ViodpllClsinV) j 

^0 / 

(OO TT 
J d\p\ J di>\pf sin^ V cos2 Qjodpl ICI sin V) + ^i2(|p| |CI sin V) 


(OO 7r 
y"d|p| Jdij\p\^ sin^i^cos^i; (^^Jo(|p||CI sinV') - |j2(|p||CI sin V) 


(OO TT \ 

sin^ tjj cos ■^i ji I CI sin 


Jd\p\ ydVlpl'sin^V- QjodpiiCi sinV) + |j2(|p||C|sinv>) 





OO TT 



+877^ I /dVlpl'sin^^ ( ijod^llClsinV') " ^J2(|p| |CI sin V') 



^0 



+327r2 Mdlpl ydVlpl^cos^V-iodpllClsinV) 

\0 

/ OO 7r 

+327r^ i Jd\p\ y#|p|^sinV'cosV'«ii(|p||CI sinV') 

D TT 

y"#|p|3 sin^ QiodpllCI sinV') + ^j2(|p||CI sinV') 



+367r2 





OO -TT 



+727rM y^lpl y#|p|3sin2^QjodPllC|sinV)-^i2dP||C|sinV) 



kO 
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We summarize the integrals of the same type and obtain Eq. (|3.30D : 



oo 



(eoo(x)eoo(y))™^ - TCT^ I 9671^ I ld\p\ IpP Id^P sin^ V-cos^ V'jodpHzl sin V) 



(27r)8/38 



OO 

+ 487r^ I Jd\p\ IppydTA sin^7Acos^'(/'j2(|p||z|sin7/;) 

\0 

oo TT \ 2 

+ 967r^ I Jd\p\ \p\^ Jdi: sin^ V cosi/'Ji (|p||z| sin V') 

\0 

^^IpI \p? jd^ sin''?Ajo(|p||z| sinV') 



oo ^ 

d\p\W jdip sin'''(/;i2(|p||z| sinV') 


oo 7r 

+ 32^^ I jd\p\\p\'^ fdijj cos' iPjo{\p\\i\sini,) 

oo 

327r'^ I Jd\p\\p\'^ J di/j sin tp cos tpji{\p\\z\ sin ^) 



oo TT 

+ 247r2 I jd\p\ \p\^ jdi^ sin^iljjoi\p\\i\ simp) 

\0 
oo 7r \ 2\ 

+ 487r2 I j d\p\\p\^ jdijj sin^ ipj2{\p\\i\simlj) 

\0 





2 



32^ 
3 

16^2 



(B.25) 



where now \p\ = \/ Pq + pf + p^ + P3 is the modulus of eucUdean four-momentum. In Eq. (j3.30p 
the last four lines arise from the term cx u^Uj/ in the propagator, see Eq. (|3.3p . Again the spherical 
Bessel functions of the first kind jojJijj2 are involved. The evaluation of integrals concerning 
functions of this kind with nontrivial argument is the subject of interest in the next subsection. 

B.3.2 Integration over first polar angle and momentum 

In order to solve the remaining integrals emerging in the vacuum contribution we use the relations 
Eqs. (|B.12p . (|B.13p and employ Sonine's first integral formula and its generalizations [60 . It holds 
that 



dV- J^(a sin V') sin^'+^ V cos'^'''^^ V = 2P+^r{p + l)a"''" Vp+^+i(a) (B.26) 
where D\efj., ^tv > — 1, a > 0, 

y"dVJM+i(asin^)sin''^ = V^2-2-Mai+Mr(i + /,),F2 J^l + 2 + J + /i; , (B.27) 
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where 9^c/x, ^Rciy > — 1 , a > 0, 



jdip sin" ^ cos" il>Jp{a sin il>) = 2-^-p{1 + {-iy)aPr 



■ 1-F2 



1 + H + p 2 + ii + v + p 



+ ) (B.28) 



provided that '^{tv > — l,9lep > — l,9le(/z + p) > —land a > 0). In these formulas pFq denotes a 
regularized hypergeometric function, which is defined as 



pFq{ai,...ap\l3i,...(3q\z) = 



pFq{ai, . . . Up] Pi, . . . (3q; z) 

r(/?i) TO 



(B.29) 



We observe that the integrals which involve the spherical Bessel functions ji (and in general 
spherical Bessel functions of uneven degree j2u-i) vanish for reasons of symmetry: 



IpI dtp sin tp cos tpj i{\p\\z\ sin ijj) = , 



oo 7r 

J d\p\ \p\^ Jdtp sin^cos V'ii(|p||z| sin^) = . 



Employing the Sonine integral formulas yields: 

oo 7r 

Jd\p\ \pf Jdtp sin^ ^cos^ Vio(|p||z| sin^) 



dip sin2 ■i/'cos^ V-^idpll^l sin^) 



oo 2 



— oo 

'^Jd\p\ |p|i(^2^r(0 (|p||z|)-ij2(|p||z|) 



I2tt 





oo 



(B.30) 
(B.31) 



(B.32) 



IpI dtp sin ■^cos Vj2(li^l 1^1 ^^'^ V') 





oo 



^j^ry Jd\p\ \p\^ jdip sin2 ^cos^'^J|(|p||z|sin^) 



= 3 





oo 



TT , _ ,_ 3 / ,1 _i I _i 3 
Z 2 



dip sin 2 ^cos xp 



|p||z| sini/) 2 



J|(|p||z|sinV') - Ji(|p||z|sinV') 



'i\p\ \p\ ^ J dip sin2 ^cos^ ■!/'J| (|p| |z| sinf/;) 



I oo f- 

2tt f 5 f 3 

— W — y d\p\ li?! 2 j dip sin 2 -^cos^ V-^idPll^l sinip) 
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CSO TT 

= Jd\p\ \p\^ JdTp sin5^cosV |^|^||^|^^T^Ji(|p||z|sinV') - J_ i (|p||z| sinV)) 





J d\p\ IpI 2 J dip sin~2 ■j^cos^ ^/'Ji (|i5||z| sintp) 



OO "2 OO 

-3\/27r|z|~^ jd\p\ \p\i J dtp sini ip cos^ 'tpj_i{\p\\z\smil)) - jd\p\ |p|J2(|p||z|) 



||z|-i /d|p||p|W2-i(i + (-i)2)rQ 



zD^ri- 



1 \ 1-^2 2' ^' 2' 



OO 

-3V2^|z|-i y"d|p| (^25r 



^I)-^^i(IpII^I) 



r(2)r(|) 



d\p\ \p\J2{\pM) 



OO 



'2z 
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-6Ji(|p||z|) - 2|p||z|J2(|p||z|) +3|p||z| 1F2 ( -; -,2; 



X 3 



00 TT 

y#i \p\' jdtp sin ^jo(|p||z| sinV)) 



CO 7r 

= y"#l \P\^ jdip smi ipj^{\p\\i\sm.ip) 



00 7r 

= y"'^|p| \P\^ jdtp sin5 V' (^j^||^^7^^|(|p||z|sinV') - J|(|p||z|sinV') 



TT 

00 2 

= 3\/27r|z|~3 y^lpl Jdtp sin^ ipJ^{\p\\z\smip) 



27r|z| 2 Jd\p\\p\^ J dip sm2 tp J ^{\p\\z\ sin tp) 


00 

= 3^|z|-i Jdlpl (2-ir (|p||z|)-i J2(|p||z|) 





' 

00 00 

^ jd\p\ \p\j2{\m\) - ^ jd\p\ \p?m\pM) 



(B.33) 



(B.34) 
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OO TT 

Jd\p\ \pf jdi, sin ^j2(|p||z| sinV') 

"2 

■/#sini^J.(|p||z|sin^) 



d\p\ \p\^ 



i\)--^M\pM) 



' 

OO 



OO TT 

jd\p\ Ipfjd^ cos VJo(|p||z| sinV)) 

fc?|p| IpI^ /df/' sin~2 ^cos^'^Ji(|p||z|sin^) 





OO 



TT 

2M 



= d\p\\p\^ 2-Hi + {-ir)r - t - 



2|z| 



1^2 (^;2, 2 



3. (IpIIzI)^ 

4 



r(2)r(|) 



c^IpI IpI /# sin -!/)jo(|p||z|sinV') 



c^IpI Jdip sin^ tpji{\p\\z\smtp) 



= f^Jd\p\ \p\i[2-iT 

' 

OO 



IpI"^ /(if/' sin ■0j2(|p||z| sin^j) 





OO 



^ yiilpl IpI^ sin^ VJidpllzlsinV-) 



D TT 

d\p\ IpI^ / df/) sin5 ip 



TT 

2|¥ 





OO 



= 3w^|z| 2 /(i|p||p|2 Idipsiniil) 



\p\\z\ sinip 
1 



27r 



|z| sin^ 



|p||z| sinf/i 



OO 2' 



J|(|p||z| sin^) — Ji (|p||z| sin^) 
Ji(|p||z|sinV)- J_i(|p||z|) 
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3y'^||z| J d\p\\p\'2 J di; sin ^Ji (|p||z| sinTA) 



oo -J 

-3^/2^|iri/.|p||p|i/d^sin^^J_.(|p||z|sin^) 



oo 

-V2^\i\-i [d\p\ \p\^ (|p||z|)-^ Ji(|p||z|) 







syflzr^ Jd\p\ ipi-i ^2-1(1 + {-in\m)ir r 

\ 

00 

-3V2^|z|-i Jd\p\ 



r(i)r(|] 



i\)-^M\p\\i\) 





00 

'•1^11^12 



z 







00 

= -Jd\p\^(3M\p\m) 



Mm\ [M\m\) - ^ (i; 1, 1 j j j (B.38) 

The analytic evaluation of the integrals over the modulus of the four-momentum is a very 
sophisticated task. In the following we compile the necessary prerequisites for a detailed calculation 
of the remaining integrals. 

At first we have to mention that the involved hypergeometric functions 1F2 (^;l,f;^), (jsfi^; 
1F2 (2;3,3;z) and 1F2 (3;4,4;z), can be expressed in terms of elementary functions and Bessel 
functions. It holds that [59] 

iJ^2 =/o(2Vi) + J(/o(2yi)Li(2V^)-/i(2Vi)Lo(2Vi)) , (B.39) 



ii^2 ( ^;|2;z') = 2/o(2VI)+7r(/o(2Vi)Li(2Vi)-/i(2VJ)Lo(2VI))-^/i(2VI), 



(B.40) 



iF2(2;3,3;z) = ^(l-/o(2^/i) + ^/I/i(2Vi)) , (B.41) 
1 8 

1F2 (2; 3, 3; z) = — ((z + 2)/o(2yi) ~ 3Vl/i(2V^) - 2) , (B.42) 

where = denotes a modified Bessel function of the first kind and represents a 
modified Struvc function. These functions are related to their regular counterparts via 

I,{iz) = i^J^iz), (B.43) 

L,(^z) = ^'^+lH,(z). (B.44) 
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The primary definition of these functions are given as power series 



The parity of Bessel and Struve functions is determined by the order v in the foUowing way: 

Iu{-z) = {-lYI,{z), (B.49) 

U{-z) = (-l)^+iL,(z), (B.50) 

Ju{-z) = {-lyj.iz), (B.51) 

H„(-z) = (-l)^+iH„(z). (B.52) 

Eqs. (??), (|B.44p and (??) - (IB.52[) show that the sign of the square root of z has no influence on 
the final result in this concrete case, since the sign cancels out anyway. Thus we constrain ourselves 
to the positive sign for the remainder of the calculation. In the following x will denote the product 
|p||z|. As a consequence of the mentioned properties of and Li,, or Ji, and H^, respectively, we 
can transform Eqs. (|B.39|1 and (|B.40ll into 

i^^sQ;!,!^) = Jo(x) + |(-Jo(x)Hi(x) + Ji(x)Ho(x)) , (B.53) 
^^'(M'^'~f") " 2Jo(x)+^(-Jo(x)Hi(x) + Ji(x)Ho(x))--Ji(x), (B.54) 



iF2('2;3,3;^] = ^ (l _ _ | , (B.55) 



-X'\ _ 64 

X' 

iF2(^3;4,4;^) = (^_^ + 2^ Jo(^) + ^ j^(^) _ 2^ . (b.56) 

Another important property of the Bessel functions is their normahzation. It holds that 

dzJ^{z) = l, (B.57) 



for positive integral and real v. 

The strategy for the evaluation of the momentum integration is to reduce every integral to a 
generalized Lipschitz-Hankel integral. An investigation concerning these integrals of the structure 

dze-"V^(/3z)z^-i (B.58) 
is presented in |601 . The interesting integrals are available by performing the limiting process 
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a ^ 0+. In this prototype integral (|B.58[) is evaluated as 



dze-'^''J^{(3z) z' 



a^'T{v + 1) 
aA'r(zy+ 1) V 

(f)'^ro. + z.) 

(a2+/32)^(A'+'-)r(i/ + l) 



2 ' 2 



i^ — /i+l V — ^ 



2F^ 



2 ' 2 
/i + i/ 1 — /i + i^ 



1: 



q;2 + 13"' 



(B.59) 



provided that 9le(^ + v) > and 9^e(^) > |Ke(i^)| to secure convergence at the origin and lHe(a ± 
i(3) > to secure convergence at infinity. The hypergeometric function 2F1 is called Gauss' 
confluent hypergeometric function. Relevant for our considerations are the cases where ^ is taken 
equal to + 1 or + 2. Then we obtain 



00 

/ 



dze-°"'J^{(3z)z'' = 
dze-°''J^{Pz)z''+^ = (2a 



(2/3rr(z. + i) 



(2/3)-r (z. + 1) 



(a2+/32) 



2 VTT 



(B.60) 



(B.61) 



so long as fHe(i^) > — 9le(z/) > —1, respectively. The two integrals above imply that ioi ^ — v 
being an even integer, an integral of the Lipschitz-Hankel type vanishes in the limit a —^ 0^ . 
At this point we have collected the tools to proceed with the explicit calculation. 

00 C30 

^ J d\p\ \p\J2m\i\) = ^,jdxx i^Mx) - Mx) 

/ \ 

00 00 

2 fdxMx)- fdxxMx) 



TT 
i?j4 



V 



2tt 

1514 



(B.62) 



2|zP 



d\p\ [ -6Ji(|p||i|) - 2\p\\i\M\m) + 3|p||i| 1F2 ( i; I 2; JJeMI 



00 00 



2ip 



dxXiF2 ( ^; ^'2;-^ 



00 



/ \ 

00 00 

2 fdxMx)- fdxxMx) 



V 



+ ^ /rfxx (2Jo(x) + TT (-Jo(x)Hi(x) + Ji(x)Ho(x)) - ^ Ji(x) 
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15^14 2|7l4 



2 jdxxMx) -TT ydxx^o(x)Hi(x)+7r j dxxMx)rio{x) - jdxMx) 



X' 

^1^ 2|z|4"^r(fc + |)r(fc + |)fc!7 '^^22^'^°^^^ 



2/0+3 



^2|z|4'' A. r (fc + 1) r (fc + 1) fc! y ^ 22'=+! -^'^^^ |Z|4 



(B.63) 



OO 

^ y^lpl (3|p||z|J2(b~||z|)-(|p||z|)V3(|p||z|)) 


oo OO 

= j|l /^XX (^-^ilx) - ^o(x)) - 1^ y"c^Xx' (^'^2(x) - ^i(x) 



^ JdxMx)--^ JdxxMx)--^ Jdxx(j^Mx)- Mx)j + Jdxx'^Mx) 



oo oo 

^4-^,JdxMx)+^, JdxxMx) = -^, 



(B.64) 



oo oo 

|| IplVsdpllzl) = ^^Jdxx' J2(X) - Ji(x) 



oo oo 



= |i /rfx Ji(x) -|i /dxxJo(x) = |i 



(B.65) 
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Jdlpl |p|3,F2 0;|2;. 



1 3 v2 





CSO 



CO 



P|4 / dxx' (2Jo(x) +7r(-Jo(x)Hi(x) + Ji(x)Ho(x)) - ^Ji(x) 



CX> OO OO CO 

2 jdxX'Mx)-^ ycixx'^o(x)Hi(x)+7r J dxx'Jiix)tioix) - ^ J dxX^'Mx) 



2 z 







=0 =0 / 



OO 



"2|z|4''£ r (fc + 1) r (fc + 1) fc! y 22'=+2 



=0 

OO 

CO f ^ \k r 2k-\-4 

-^^T^y^ dx^^7TT-TJi{x)=^ (B.66) 



CO OO 

|| y d\p\ |p|Vi(|p||z|) = ^ jdxx'H\m\) = (B.67) 







^ Idlpl (^3|p|Jo(|p||z|) + |pn^|Ji(|p||z|)-|pUF2Q;l,|-^^^) 

OO 

^r(i)r(|) (2;!' 2' 4 



^ J dxxMx)--^ J dxX^Mx)+-^ ^^^^^^ss J dxXiF2l-,l,-,-^ 



=0 =0 



= J^lJdXX (Mx) + ^(-Jo(x)Hi(x) + Ji(x)Ho(x))) 



00 00 00 

= 1^ jdxxMx)--^^ jdxxMx)'H.i{x) + j dxxJi{x)^o{x) 



=0 



=0 

TT (-1)'= 7 y2'=+2 
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Finally we calculate the order-of-magnitude estimate for the vacuum contribution to the two- 
point energy density correlator in a more detailed manner. Recall that the integrations are per- 
formed up to a limit determined by the maximal resolution of the effective theory. In principle, 
we obtain our results by employing the following master integrals (with the integration cut-off A) 







where fHe(a + [3) > —I and > 0, 

dpp'^,F2ia;b,c;ap^) = ' ^ ' ' ' ^ ' ^, B.70 

J i + a 



where [He(a) > — 1 and crA^ G IR<o- An important formula concerning integrals over a separate 
Bessel function is given by the sum 

/dz J,(z) = 2^ J,+2fc+i(A), (B.71) 
fe=o 

provided that 9^e(i^) > —1. The regular Bessel functions of the first kind of even order are 
resummable into the expression 

oo ^ 

Y.J2k{z)^-{l^Mz)). (B.72) 

k=l 

In several cases we will be confronted with integrals containing a product of Bessel and Struve 
functions. Then the following indefinite integral is employed j59j : 

dz ^ = — ■ — - J^(z)H^_i(z)z 

lF2{ -,1^+1 



0Fr(Ai + i)r(;. + i)r(ii±p) ' ' V 2 ' ' ' 2 '4 

-(zJ,_i(z) + iy)J,iz))B.^{z)^ . (B.73) 

Finally we should mention the primary definition of a generalized hypergeometric function. 
pFq{ai, . . . , flp; 61, ... , bq] z) is defined in terms of a power series: 

pFq[ai,...,ap;bi,...,bq;z) = 2^ \ T\ ' ^^■^'^> 

where {a)k denotes the Pochhammer symbol 

T{a + k) 
r{a) 

This definition implies the following property: 

pFq{ai, . . . , ftp; 61, . . . , bq^i,a.p\ z) Fg_i(ai, . . . , ap_i; 61, ... , 6g_i; z) . (B.76) 

With these formulas we are able to proceed with a detailed calculation of the integrals. 
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^ Jd\p\\p\M\pM 



TT \4>\*\Z 





1-F2 2; 3, 3; — - 



1 - Io{l\(j)\\z\) + — — /i(i|0||z|) 



|Z|2 32 |<^|4|5 

= ^ (2 - 2Jo(|^||z|) - |^||z| Ji(|^||z|)) (B.77) 



101 

^dlpl f-6Ji(|p||z|)-2|p||z|J2(|p||z|)+ ii^2fi;^,2: 





I0l|i| 101 101 

/■...w.A ^ /'.I.M.I r /,.n.n , 377 Z-^,^, ^ /I 3 ^_ |p|2|z|^ 



/ dx Mx) - ^ /rflpl \p\J2{\pM + ^ jd\p\ 1F2 Q; ^, 2; - 
00 

l^zf; J2(,+i)(x) + ^ (2Jo(|^||z|) + |^||z| Ji(|.^||z|) - 2) 



101 





37r 



(1 - Jo(|0||z|)) + ^ (2Jo(|.^||z|) + |^||z| Ji(|.^||z|) - 2) 



.^HW..3 ii;i,^,2;-Mi __|^(,_.,,^,,,)) 



^ ( 2 (8Jo(|0||z|) + |<^||z| Ji(|<^||z|) - 8) + 3|^nzp 1F2 ( 1; | 2; 



101 



(B.78) 



z 





101 101 

^ j d\p\ \p\J2{\pM) - ^ / d\p\ |p|V3(|p||z|) 



Mmz\) - -^JmM) j - (^3;4,4; — j 



|z| 288 (-|^|6|z|6) 4 J 2 

-W (^'*^'''^'' ~ 2) Jo(|.^||z|) + 3|.^||z| Ji(|^||z|) - 2) 



(B.79) 
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101 

d\p\ WM\pM) 







1^ ((|^nz|2 - 2)Jo(|^||z|) - 6|^||z| Ji(|.^||z|) + 2) (B.80) 



101 



4 





2 4 M 2' ' 2' ' ' 4 



101 





1011^1 

2 



|z|4 / 


= ^(|^nz|Vi(|^||z|)) (B.82) 



101 

^ y^lpl (3|p|Jo(|p||z|) + |pn^|Ji(|p||z|)-3|pUF2Q;|2;. 



|p|2|z| 







1011^1 1011^1 101 

y c^XXJo(x)-p|4 y rfxx'^i(x)+ ||p J d\p\^F2i^-;-,2;- ' ' 
000 



|Z|4 



-|il<^||z| Ji(H|z|) - ^l^l^lzl V,(|^||z|) + ^ ,F3 (^1, 1; 1, 1 2; j 

|^6|^||z| Ji(l^lizl) + 2\4>\'\ifj2{\m) - m'l^f 1F2 (^l; I 2; j j 

(B.83) 
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B.4 Calculation of the vacuum and thermal contribution in 
Lorenz gauge 

We evaluate the two-point correlator of energy density again in Lorenz gauge in order to show 
the gauge invariance (in case of the vacuum contribution) of this quantity explicitly. As men- 
tioned above, we have used Coulomb gauge {Qoo{x)Qoo{y)) in the case of the 5[/(2)-Yang-Mills 
thermodynamics calculation, since this gauge is a physical gauge and the implementation of the 
constraints on the momenta is given as in Eqs,. (|2.96p . (|2.97p . and thus considerably easier than in 
the covariant gauge. 



B.4.1 Propagator of the massless gauge mode in different gauges 

For reasons of completeness we will compile the the propagators of the massless gauge mode 
in momentum space and position space in consideration of different gauge conditions. We will 
constrain ourselves to the cases of Lorenz and Coulomb gauge. We recall that in the real-time 
formalism of finite temperature field theory the propagator consists of a vacuum (zero-temperature) 
and a thermal contribution: 

D^Ap) = D-{P, T) + D%{p, T). (B.84) 

For a more detailed discussion one could consult for instance [ID]. 
Coulomb gauge: 

The momentum-space vacuum contribution is given as 

D-[P) = -Pl-^ + - , (B.85) 

p + le p Pq 

with the gauge-fixing parameter. The quantity {Qooix)&ooiy)) does not depend on any value of 
the gauge parameter ^. The usual Coulomb gauge is realized for the case ^ = 0. The contribution 
coupling to the zero-component of the electromagnetic current is the Fourier-transform of the 
Coulomb-potential. In position-space we obtain 

D;7-~{\^ y|) = -5{xo - yo)^^ , (B.86) 

and describes an instantaneous static interaction. The momentum-space thermal contribution is 
given as 

p) = -pj;(p)27r^(/)nB(/3boi) . (b.s?) 

Lorenz gauge: The momentum-space propagator for the massless gauge mode in covariant gauge 
is given as 

D^Ap) = -9,^^h- + ^ 7 \ ^2 - 9M'^^Kp>b{(3\po\) ■ (B.88) 

p + £, (P + *£) 

Here ^ again denotes the so-called gauge-fixing parameter. A convenient choice is — 1, which 
is referred to as Feynman gauge, or the Landau gauge ^ = oo. We calculate in Feynman gauge. 
While apart from the tensor structure the thermal part is equal to the Coulomb gauge thermal 
part, the Lorenz gauge vacuum propagator in position gauge is then given as 

D-ix -y)^ {T[A,{x)AAv)]) = • (B.89) 

In one of the following subsections we will use this propagator Eq. (jB.89p in order to proove the 
above obtained result for the vacuum part by a straight-forward position-space calculation. 

Via employing a Fourier transform we are able to determine position-space counterpart of the 
thermal propagator. We obtain the equal-time correlator 

/3-^(x-y)coth(i(^) 
^-(^ - ^' ^) ^..0 = -.9- 2./3(x-y)^ • (^-^^^ 
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A more sophisticated task is the calculation of the nonequal-time thermal part of the propagator. 
Nevertheless by consulting [H] we are enabled to evaluate the four-dimensional Fourier integral, 
which yields 




sinh 



27r|x-y| 



/3|x - y| cosh (^^^f^) - /?|x - y| cosh (^^iIblzM^ 7r|x - y|2 - 7r|a;o - yo| 

(B.91) 

In the hmit x° y° Eq. ((R9T|) is equal to Eq. ((R90l) . 



B.4.2 The thermal contribution: Calculation in position space 

The thermal part can also be computed in coordinate space. Instead of calculating integrals, we 
insert the propagator Eq. (|B.9ip into Eq. (|B.2p and calculate derivatives of this structure. Finally, 
we set = = and obtain the equal-time two-point correlation of energy density. 
Nevertheless, the derivatives of Eq. (jB.9ip with respect to x and y are very large expressions and 
involve hyperbolic functions. Although these derivatives can be simplified by making use of several 
functional identities, a calculation in momentum space is considerably more convenient. 
In order to compile the necessary ingredients for a momentum space calculation of (6oo(2;)6oo(y))"' 
we give the following formulas: 



coth (^) cosech^ (iM) 

/ 



67r2|zp/?- 37r|z|/32 sinh 



27r|z| 



ZiZn cosech^ 



^^j^^oc...! I /3 y /,_3 . 2 / ttIzi 



. , m^. I le^sinh" (B.93) 

9,o5,,i?-(0,/?)|^„^g = d,,dyoD'^,{z,p)\^,^^ = 0, (B.94) 
where z stands for x — y. 

B.4.3 The thermal contribution: Calculation in momentum space 

Let us first consider the thermal part to the correlator. We insert the propagator 

D^^^{p,T) = ~g^,{p)2Ti5{p^)nB{P\po\) (B.95) 

into the Wick decomposition of the energy density two-point correlation function, Eq. (|B.2p . Sub- 
sequently, we rearrange terms of the same structure and use the properties of the Delta distribution 
Eqs. ()3.11l [3TT2|) and integrate over the zero-coordinates. This yields 
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(eoo(cc)eoo(y)>"' 

•(5(po - IpI) + Sipo + IplWiko - |k|) + Siko + |k|)) 



Q-ipix-y) ^-ikix-y) 



iS{po - IpI) + S{po + |p|))(<5(fco - |k|) + S{ko + |k|)) 

//" f d?P f d?k 

■iSip„ - IpI) + ,5(po + |p|))(^(fco - |k|) + 5{ko + |k|)) 

•(<5(po - IpI) + 5{po + |p|))Wfco - |k|) + 5{ko + |k|)) 



<Pp f <Pk /„,.x2.„^//,i IN ^//oi, |N.^-ip(x-y)^-ife(x-y) 



+ '^jdpoj dko j ^2Trf 2\p\ J (27r)32|k| ^^^^ ' "^(/^bol) "s(/3|A;o|) • e-'J'^^-^^e 
•(^(Po - IpI) + S{po + |p|))(<5(fco - |k|) + S{ko + |k|)) 

- 2^00 / dpo I dko J j^^^ I J^^^iPkr ■ nsiPlpol) nB(/3|/co|) • e--(-^)e-'=(-^) 

■{S{po - IpI) + S{po + |p|))(<5(fco - |k|) + 5{ko + |k|)) 

- 2ffoo J dpo J dko J J (2j^)32|k| ^°^' • "^('^l^oD ^B(/3|fco|) • e-^^'(--'')e-^'=(--'') 

•(<5(Po - |P|) + S{po + \p\)){6{ko - |k|) + Siko + jkj)) 

- 2,oo/cipo/*o/(^/(^/fcg-nB(/3bo|)^ 

•('5(po - IpI) + S{po + |p|))(<5(fco - |k|) + <5(fco + |k|)) 
+ glo j dpo j dko j (^Jj32|p| / (2j^32|k| (P^)' • "^(^l^'ol) «B(/3|fco|) • e-^'^^-^'e-'^^--^) 

•('^bo - IpI) + 5{Po + |p|))Wfco - |k|) + 5{ko + |k|)) 
+ f ^^o/.fco/^/(^/.^-n.(/9|po|)n.(/^|.o|)-^^^^^^^^^ 

•(^(Po - |p|) + S{po + \p\)mko - |k|) + 5{ko + |k|)) 

(B.96) 
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•(cos((|p| + |k|)(a; - y)o) + cos((|p| - \k\){x - y)o)) 
•(cos((|p| + |k|)(:r - y)o) + cos((|p| - \k\){x - y)o)) 

d^p f d^k ^ Pk ^ |p||i^|.^^(^|p|)^^(^|i^|).^,p(x-y)^,k(x-y) 

(cos((|p| - mx - y)o) - cos((|p| + |k|)(x - y)o)) 



27r)3 J (27r)3 V|p||k| 



+1 



■(co8((|p| + |k|)(x - !,)„) + co8((|p| - IklKi - y)„)) 

d^p f d^k f pk 



(2^)3 y (2^)3 v|p||k| 



, d^p f d^k ( 

-900 / 7K-^ / 7K-^ 1 



(27r)3 J (27r) 



p||k| • nsim) nsim) ■ e'P^""''' e^''^^-^' 
•(cos((|p| - |k|)(:r - y)o) - cos((|p| + |k|)(x - y)o)) 

\p\\k\-nBmp\)nBmk\)-e^P<^--yU^^(-y^ 
•(cos((|p| + |k|)(:r - y)o) + cos((|p| - |k|)(x - y)o)) 
IWflWf IPllkl •-B(/3|p|)n.(;3|k|) .e-(>^-)e-(>^-) 

•(cos((|p| - \k\){x - y)o) ~ cos((|p| + |k|)(x - y)^)) 

•(cos((|p| + |k|)(a: - y)o) + cos((|p| - \k\){x - y)o)) 
(cos((|p| - \k\){x - y)o) - cos((|p| + \k\)ix - y)o)) 



-5oo 



5oo 



/■ d3fc 



(27r)3 J (2 



^ (l+ (^)') |p||k| •nB(/3|p|)nB(/3|k|) .e*P(--^)e''^(''-''' 



(B.97) 



Last line is equal to the expression derived in Coulomb gauge, Eq. p.l3p . Subsequently we can 
proceed in exactly the same manner as we have done above. 

B.4.4 The vacuum contribution: Calculation in position space 

If we would only were interested in the value of (6oo(2;)6oo(y))™'' a calculation in position space 
and Lorenz gauge would have been more convenient. For the reasons mentioned above we have 
preferred a calculation in momentum-space, in order to present how our method works, particularly 
with regard to the more involved calculation, when photon propagation is governed by an SU{2) 
gauge principle. At this point we briefly compile the fundamental elements of the position-space 
computation of the energy density two-point-correlator and surrender details, since the calculation 
is of straight-forward character and more a challenge on concentration. 

Starting point for the calculation is again the Wick decomposition of the correlator, Eq . (|B.2[ ). 
Instead of shifting to momentum-space, we directly insert the position-space propagator Eq . (|B.89p . 
Differentiation with respect to the x- and y-coordinates and contraction of the tensor structure 
will yield us the final result. 
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Concerning the calculation we have to consider derivatives of the propagator. It holds that 
{T[d,.A^{x)d,.Ax{v)]) - d,.d,.{T[A,{x)Ax{y)]) 

= 5«A-^(^ j,(B.98) 

where T denotes the time-ordering symbol. Since negligence of time-ordering is well justified we 
omit the symbol in the following. In principle, we have to evaluate two-point correlators of squared 
field strength tensors. It is not difficult to see that following relation holds: 

{T[F\x)G\y)]) = 2 {T[F{x)G{y)]) (B.99) 

The correlator of field strengths is given as: 

{T{Ff,^{x)F^x{y))) = 9t,iJ^u\{x -y)~ g^\I,,n{x - y) - g,yf,I^\{x - y) + g,yxlf,^{x -y){x-y). 

(B.lOO) 

The energy density correlator is then given as the sum 

(eoo(a;)eoo(2/))™^ = 2To^ {x - y)Toxor{x - y) - ^""(a^ - y)J'oxar{x - y) 

(B.lOl) 

Inserting Eqs. (jB.QSp and (jB.lOOp and contracting contravariant and covariant tensor components 
let us finally arrive at 

(eoo(a;)eoo(2/)) = tttTZ — 7\K 3.9oo - 8goo^ 7^ + 32 



27r4(a; - 


-yf * 


1 




47r4(x ^ 


„y)8 ( 


1 




47r4(x ^ 


-yf ' 


1 




7r4(a: - 


yf ( 


1 


( 



{x - y)2 {x - y) 

2 ic {xo-yaf 



[x - yy 

2 /in (^0 " 2^0) 



-11.900 - 40goo-7 

2 (xq — yo)^ 

22^00 - IO500 



[x - yy 

Assuming equal time and setting x° — we obtain exactly the same result as in the previous 
calculation in Coulomb gauge. 

B.4.5 The vacuum contribution: Calculation in momentum space 



Again Eq. (|B.2p provides a starting point for the calculation. Insertion of the propagator in covari- 
ant gauge yields 

(eoo(a;)eoo(2/))"' = 

' (4p2fc2 - Apokop^k" + 2plk^k'' + 2klp^p^' + 2p^p^'k,k'') 



{2ttY J (2^)4 ^ " ^ '^"'^^ ' p^+ is fc2 + j£ 

l^'P I {p,p^k,k^ + Aplk,k^ + Aklp,p^)'' 



4 J (27r)4 J (27r)4 ^^^^ ^" "^^^ ^ p2 + ie fc2 + 

goo I dp I dk ^p^p.k^k'^^iplk.k^ + ik^p.p'^)'' 



4 J (27r)4 7 (2^)4 ^^^^ ' 'J^'^^ ^ p2 ^_ j£ fc2 + 

o2 /■ rf4 f Jih. p-ip{x-y) p~ik(x-y) 

9oo I I (16p^p^fc,fc^ + Sp^p^^k^kn % , . ^T^— . (B.103) 



16 y {27r)4 J (27r)4 ^ " > p2 + + ie 
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By virtue of analytic continuation we rotate to euclidean signature and arrive at 
(eoo(a;)eoo(2/))* = 



(27r)4 J (27r)4 ^i^u.^u^^.^^ ' -^u^/x.^^ , -.^u^^m^m ■ -.-m^m— y ^2 ^2 

- J y {Pi.Pi.k.k^ + "^PoKK + 4fcoW) ^2 p— 

+ J J {16p^,p^k^k^ + Sp^p^k^k^) —^^ ^ . (B.104) 

Finally we obtain 

d^p f d^k ( \\ e*P(="-2') e^^'^-v) 



p||k|y J p2 1^2 



+ 5 TT^^Po 



d^p 2 e^P'^^'-y^ 



2 



• fi^p /• / ^ pk gip(x-y) gifc(x-j/) 



p 



n2 



2 



2 Vy (27r)4 p' 
/ /■ d*n c'^Pix-y)\ / f d'^k Qik{x-y) 



rf^P , 2 e'P(^-^') \ ( [ d'^k 2 e'*=(^-2') ' 

(B.105) 



In the following we proceed exactly in the same manner as in the case of the Coulomb gauge 
calculation. For simplicity we set x° = y° = 0, rescale to dimensionless variables, introduce hyper- 
spherical coordinates and perform the integrals successively. 
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(eoo(a;)eoo(2/)) 
1 



(27r)8/38 



6 • 27r y dlpldippclBj, J fi|fc|rf^ferf^fe|p|^ sin^ cos^'p sin^p|fc|^ sin^ Vfc cos^'fe sin^fe 

_ gi|p||z| sini/'p cosepgi|fe||z| sini/'fe cosOk 

■ / (i(^p (sin 9p sin (^fc cos (pp + cos 0p cos 9k) 
Jo 

+5 • 27r / (i|p|d'!/'pC?^p / (i|A|(i'i/'fe(i6'fc|p|'^sin'*?/'pSin(?p|fc|^sin'''!/;/5sin0fc 



. g»|p|l^l sini/'p cos6pgi|fc||z| sini/'fe cosSfe 
r2v 



■ I d(f)p (sin^ 9p sin^ 9k cos^ (j)p2 sin sin 9k cos cos cos (j)p + cos^ cos^ 9k) 
Jo 

+ i -477^ (i|i>|d'!/'pd0p|p|^sin^'!/jpCOS^'!/'pSin0pe*' 

• d\k\dijjkd9k\k\^ sin"^ tpk sin 9 ke''^ 

Y d|fc|#fcd6'fe|fc|3sin2V'fcCos2V'fcsin6'fee'l^ll^l'^^'^'^'='=°"^'=^^ (B.106) 



^i|/c||z| sintpk cosOk 



OO TT +1 

G-Att^ ljd\p\ Jdip y"dcos6'|p|3sin3V'cosV'Pi(cos6')e*IPll^l^'"'''^°^^ 



0-1 

2 



(27r)8/3' 

(OO TT +1 

y"d|p| y"d?A Jdcos9\p\^sm'^ipcos'^ 'iIjPo{cos9) e'l^l 1^1 ^i"V' cose 
0-1 

(OO TT +1 

Jd\p\ Jd^J JdcOs9\p\^Sm'^i;(^^Po{cOS0)-'^P2{cOs9)^ gi|p||i|sini/,cos9 
0-1 

fOO TT +1 ^2 

Jd\p\ Jdip Jdcos9\pfsm'^i;(^^Po{cos9) + ^P2{cos9)^ ^i\msini>cos9 



1 9 

._ . 47r2 



VO -1 

OO TT +1 



Jd\p\ JdiP y"dcos0|p|3sin^V^b(cos^)e*l^*ll^l"'"^'=°^^ 

0-1 

OO TT +1 

+4-47r^ iyd|p| J dipp J d cos 9p\p\^ siio^ tppcos^ tppPo{cos9p) 

0-1 

OO -TT +1 

d|fc| jdi^k jdcos9k \k\^ sir? Vfefo(cos6'fe) e'l^H^I 



i|p| |z| sin i/jp cos 6p 



\0 
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CSO TT +1 



KO -1 

oo IT +1 



<0 -1 



/ 



(27r)8/38 



oo TT 



487r I J d\p\ yd'i/'lpl sin '(/'COS'(/'iji(|p||z|sin'!/;) 





V 



=0 



oo TT 

+807r^ ( Jd\p\ JdifjlPl^ sm^ ip cos^ tpjo{\p\\z\smip) 

\0 

y"d|p| y"#|p|^sin^V (^^Jo(|p||z|sinV') + |j2(|p||z|sinV') 



VO 

CO TT 



+807r^ I ld\p\ /dVlpl^sin^V ( ijo(|p| |z| sin V') - ^j2(|p||z| sinV') 



VO 

oo 7r 



sin Vjo(|p||z| sinV)) 



\0 



=0 



oo 7r 



oo 7T 



+647r /dVpIPl sin ^pcos Vpjo(|p||z| sin'0) sin Vfcjo(|A;| |z| sin V'fc) 



(B.108) 



vo 







oo TT 



oo TT 



+ 647r^ 



pIpI^ sin^ ^/ipjodpllzl sinf/') j ( Jd\k\ Jdipk\p\^ siu^ ipkcos^ tpkjo{\k\\z\smtpk) 



^0 







(27r)8/38 



=0 

oo TT 



807r^ 



ld->p\p\^ sin2 i/jcos^ V-^i (|p||z| sin 







oo 7r 



leoTT^ / 1 



J dtp \p\^ sin ^ il; J i{\p\\z\ sin ip) 





oo TT 



+ - 



leoTT^ / 1 



/2|z| 



d\p\ J dip \p\^ sin ^ ip J i{\p\\z\ sin ip) 





+ 



807r2 / 1 



OC' TV 



+- 



9 

3207r2 / 1 



Jdtp \p\ 2 sin 2 ^Ji (|p||z| sin'^) 





oo TT 



d\p\ I dip \p\^ sin 2 %pj 5 {\p\\z\ sin ip) 
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(B.109) 



Exploiting the integral formulas which we have proved in the previous Sec. (|B.3.2p yields: 

2tiV leOTT^ /-27r' " 



(27r8)^8 I \^|2|4 J 9 

leOTT^ / Stt \ ^ SOtt^ / -27r\ ^ 32071^ / Stt ^ 



+- 

15 



9 



(B.llO) 



This results is exactly equal to the result obtained in Coulomb gauge. This is another evidence for 
the gauge invariance of the two-point correlator of energy density. 
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